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Abstract. Luminous BA-type supergiants have enormous potential for modern astrophysics. They allow topics ranging from 
non-LTE physics and the evolution of massive stars to the chemical evolution of galaxies and cosmology to be addressed. A hy- 
brid non-LTE technique for the quantitative spectroscopy of these stars is discussed. Thorough tests and first applications of the 
spectrum synthesis method are presented for the bright Galactic objects 77 Leo (AO lb), HD 1 1 1613 (A2 Iabe), HD 92207 (AOIae) 
and /30ri (B8 Iae), based on high-resolution and high-S/N Echelle spectra. Stellar parameters are derived from spectroscopic 
indicators, consistently from multiple non-LTE ionization equilibria and Stark-broadened hydrogen line profiles, and they are 
verified by spectrophotometry. The internal accuracy of the method allows the lcr-uncertainties to be reduced to <l-2% in T e g 
and to 0.05-0.10 dex in logg. Elemental abundances are determined for over 20 chemical species, with many of the astrophysi- 
cally most interesting in non-LTE (H, He, C, N, O, Mg, S, Ti, Fe). The non-LTE computations reduce random errors and remove 
systematic trends in the analysis. Inappropriate LTE analyses tend to systematically underestimate iron group abundances and 
overestimate the light and a-process element abundances by up to factors of two to three on the mean. This is because of the 
different responses of these species to radiative and collisional processes in the microscopic picture, which is explained by fun- 
damental differences of their detailed atomic structure, and not taken into account in LTE. Contrary to common assumptions, 
significant non-LTE abundance corrections of ~0.3 dex can be found even for the weakest lines (W A < 10 mA). Non-LTE abun- 
dance uncertainties amount to typically 0.05-0. 10 dex (random) and ~0.10dex (systematic lcr-errors). Near-solar abundances 
are derived for the heavier elements in the sample stars, and patterns indicative of mixing with nuclear-processed matter for the 
light elements. These imply a blue-loop scenario for 77 Leo because of first dredge-up abundance ratios, while the other three 
objects appear to have evolved directly from the main sequence. In the most ambitious computations several ten-thousand spec- 
tral lines are accounted for in the spectrum synthesis, permitting the accurate reproduction of the entire observed spectra from 
the visual to near-IR. This prerequisite for the quantitative interpretation of intermediate-resolution spectra opens up BA-type 
supergiants as versatile tools for extragalactic stellar astronomy beyond the Local Group. The technique presented here is also 
well suited to improve quantitative analyses of less extreme stars of similar spectral types. 
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1. Introduction 

Massive supergiants of late B and early A-type (BA-type su- 
pergiants, BA-SGs) are among the visually brightest normal 
stars in spiral and irregular galaxies. At absolute visual magni- 
tudes up to My ^ -9.5 they can rival with globular clusters and 
even dwarf spheroidal galaxies in integrated light. This makes 
them primary candidates for quantitative spectroscopy at large 
distances and makes them the centre of interest in the era of ex- 
tragalactic stellar astronomy. The present generation of 8-10m 
telescopes and efficient multi-object spectrographs can poten- 
tially observe individual stars in systems out to distances of the 

Send offprint requests to: N. Przybilla (przybilla@sternwarte.uni- 
erlangen.de) 

* Based on observations collected at the European Southern Ob- 
servatory, Chile (ESO N° 62.H-0176). Table lA~T1 is available in elec- 
tronic form only. 



Virgo and Fornax cluster of galaxies (Kudritzki 1998 2000 1. 
The first steps far beyond the Local Group have already been 
taken (Bresolin et al. l200T1l2002l Przvbilla l2002l i. 

BA-type supergiants pose a considerable challenge for 
quantitative spectroscopy because of their complex atmo- 
spheric physics. The large energy and momentum density of 
the radiation field, in combination with an extended and tenu- 
ous atmosphere, gives rise to departures from local thermody- 
namic equilibrium (non-LTE), and to stellar winds. Naturally, 
this makes BA-SGs interesting in terms of stellar atmosphere 
modelling and non-LTE physics, but there is more to gain from 
their study. They can be used as tracers for elemental abun- 
dances, as their line spectra exhibit a wide variety of chemi- 
cal species, ranging from the light elements to a-process, iron 
group and s-process elements. These include, but also extend 
the species traced by Hn-regions and thus they can be used 
to investigate abundance patterns and gradients in other galax- 
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ies to a far greater extent than from the study of gaseous neb- 
ulae alone. In fact, stellar indicators turn out to be highly 
useful for independently constraining (Urbaneja et al. [2005 1 
the recently identified systematic error budget of strong-line 
analyses of extragalactic metal-rich H n regions (Kennicutt et 
al. 120031 Garnett et al. 120041 Bresolin et al. [20041 and its 
impact on models of galactochemical evolution. BA-type su- 
pergiants in other galaxies allow the metallicity-dependence of 
stellar winds and stellar evolution to be studied. In particular, 
fundamental stellar parameters and light element abundances 
(He, CNO) help to test the most recent generation of evolution 
models of rotating stars with mass loss (Heger & Langer 2000 
Meynet & Maeder l2000ll2003ll2005l Maede r & Me ynet l200ll 
and in addition magnetic fields (Heger et al. 2005, Maeder & 
Meynet 2005 1. These make predictions about the mixing of the 
stellar surface layers with nuclear processed matter which can 
be verified observationally. Moreover, BA-SGs can act as pri- 
mary indicators for the cosmological distance scale by appli- 
cation of the wind momentum-luminosity relationship (WLR, 
Puis et al l 19961 Kudritzki et al.[T999 ) and by the flux-weighted 
gravity-luminosity relationship (FGLR, Kudritzki et al. 2003 
Kudritzki & Przybilla 2003 ). In addition to the stellar metallic- 
ity, interstellar reddening can also be accurately determined, so 
that BA-SGs provide significant advantages compared to clas- 
sical distance indicators such as Cepheids and RR Lyrae. 

Despite this immense potential, quantitative analyses of 
BA-SGs are scarce. Only a few single objects were studied 
in an early phase; several bright Galactic supergiants, prefer- 
entially among these aCyg (A2Iae) and 77 Leo (Groth fl961l 
Przybylski [T^69l Wolf WH\ Aydin \WT% and the visually 
brightest stars in the Magellanic Clouds (Przybylski 1968 
[T97T1[T972l Wolf fT972lfT973l . Near-solar abundances were 
found in almost all cases. Surveys for the most luminous 
stars in the Local Group galaxies followed (Humphreys 1980 
and references therein), but were not accompanied by detailed 
quantitative analyses. These first quantitative studies were out- 
standing for their time, but from the present point of view 
they were also restricted in accuracy by oversimplified anal- 
ysis techniques/model atmospheres, inaccurate atomic data 
and the lower quality of the observational material (photo- 
graphic plates). Non-LTE effects were completely ignored at 
that time, as appropriate models were just being developed (e.g. 
Mihalas 1978 and references therein; Kudritzki l 19731 . 

BA-type supergiants have become an active field of re- 
search again, following the progress made in model atmosphere 
techniques and detector technology (CCDs), and in particular 
the advent of 8-10m-class telescopes. In a pioneering study 
by Venn dl995al Q995b ) over twenty (less-luminous) Galactic 
A-type supergiants were systematically analysed for chemical 
abundances, using modern LTE model atmosphere techniques, 
and non-LTE refinements in a few cases. These indicated near- 
solar abundances for the heavier elements and partial mixing 
with CN-cycled gas. A conflict with stellar evolution predic- 
tions was noted, as the observed high N/C ratios were realised 
through carbon depletion and not via the predicted nitrogen en- 
richment. Later, this conflict was largely resolved by Venn & 
Przybilla (2003). However, an analysis of helium was not con- 
ducted, and more luminous supergiants, which are of special 



interest for extragalactic studies, were missing in the sample. 
Similar applications followed on objects in other spiral and 
dwarf irregular galaxies (dlrrs) of the Local Group (McCarthy 
et al. [T995l Venn fT^l Venn et al. I2UU01 l2"0Tm l2U0lt and the 
nearby Antlia-Sextans Group dlrr galaxy Sextans A (Kaufer 
et al. 2004 1. The primary aim was to obtain first measurments 
of heavy element abundances in these galaxies. Good agree- 
ment between stellar oxygen abundances and literature values 
for nebular abundances was found, with the exception of the 
dlrr WLM (Venn et al. 2003 1. The number of objects analysed 
is small, thus prohibiting the derivation of statistically signifi- 
cant further conclusions. Moreover, abundances of the mixing 
indicators He, C and N were not determined in almost all cases. 

Parallel to this, a sample of Galactic BA-SGs were stud- 
ied by Takeda & Takada-Hidai (2000 and references therein) 
for non-LTE effects on the light element abundance analy- 
ses which confirm the effects of mixing in the course of stel- 
lar evolution. However, a conflict with the predictions of stel- 
lar evolution models was found, because C depletion was ap- 
parently accompanied by He depletion. The authors regarded 
this trend not being real, but indicated potential problems with 
the He 1 line-formation. On the other hand, stellar parameters 
were not independently determined in these studies, but esti- 
mated, which can result in severe systematic global errors as 
will be shown later. Verdugo et al. dl999b> concentrated on 
deriving basic stellar parameters for over 30 Galactic A-type 
supergiants, assuming the validity of LTE. Independently, the 
Galactic high-luminosity benchmark ffCyg was investigated 
for elemental abundances by Takeda et al. ( 1996 1, solving the 
restricted non-LTE problem, and by Albayrak (2000), using a 
pure LTE approach. Fundamental parameters of aCyg were 
determined by Takeda ( 1 994 ) and by Aufdenberg et al . ( 2002 1 . 
At lower luminosity, 77 Leo was the subject of a couple of such 
studies (Lambert et al. 119881 Lobel et al. 11992b . The results 
from these studies mutually agree only if rather generous er- 
ror margins are allowed for, implying considerable systematic 
uncertainties in the analysis methods. Among the late B-type 
supergiants the investigations focused on /?Ori (Takeda 1994, 
Israelian et al. [T997). A larger sample of Galactic B-type su- 
pergiants, among those a number of late B-types, was stud- 
ied by McErlean et al. (1999) for basic stellar parameters 
and estimating elemental abundances, on the basis of unblan- 
keted non-LTE model atmospheres and non-LTE line forma- 
tion. The chemical analysis recovers values being consistent 
with present-day abundances from unevolved Galactic B-stars, 
except for the light elements, which again indicate mixing of 
the surface layers with material from the stellar core. 

The basic stellar parameters and abundances from modern 
studies of individual objects can still be discrepant by up to 
~20% in r e ff, ~0.5dex in logg and ~1 dex in the abundances. 
This is most likely to be the result of an interplay of system- 
atic errors and inconsistencies in the analysis procedure, as we 
will discuss in the following. In order to improve on the present 
status, we introduce a spectrum synthesis approach for quanti- 
tative analyses of the photospheric spectra of BA-SGs. Starting 
with an overview of the observations of our sample stars and 
the data reduction, we continue with a thorough investigation 
of the suitability of various present-day model atmospheres 
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Fig. 1. Spectra of the sample stars around Hy and Ha. The major spectral features are identified, short vertical marks indicate 
Fen lines. Luminosity increases from bottom to top. Likewise, Ha develops from a pure absorption line into a P-Cygni-profile 
because of the strengthening stellar wind. Note the marked incoherent electron scattering wings of Ha in the most luminous 
object. 



for such analyses, and their limitations. General aspects of our 
line-formation computations are addressed in Sect.|4l while the 
details are summarised in AppendixlAl Our approach to stellar 
parameter and abundance determination is examined and tested 
on the sample stars in Sects. [5]and[6] The consequences of these 
for the evolutionary status of our BA-SG sample are briefly 
discussed in Sect. [7] Finally, the applicability of our technique 
is tested for intermediate spectral resolution in Sect. [8] At all 
stages we put special emphasis on identifying and eliminating 
sources of systematic error, which allows us to constrain all 
relevant parameters with unprecedented accuracy. We conclude 
this work with a summary of the main results. 

We will address the topic of stellar winds in our sample of 
BA-SGs separately, completing the discussion on the analysis 
inventory for this class of stars. Applications to Galactic and 
extragalactic BA-SGs in Local Group systems and beyond will 
follow. Note that the technique presented and tested here in the 
most extreme conditions is also well suited to improve quanti- 
tative analyses of less luminous stars of similar spectral classes. 



2. Observations and data reduction 

We test our analysis technique on a few bright Galactic su- 
pergiants that roughly sample the parameter space in effective 
temperature and surface gravity covered by future applications. 
We chose the two MK standards 77 Leo (HD 87737) and /30ri 
(HD 34085), the brightest member of the rich southern cluster 
NGC 4755, HD 1 1 1613, and one of the most luminous Galactic 
A-type supergiants known to date, HD 92207, for this objective. 



For 77 Leo, HD 111613 and HD 92207, Echelle spectra us- 
ing FEROS (Kaufer et al. fT999b at the ESO 1.52m telescope 
in La Silla were obtained on January, 21 and 23, 1999. Nearly 
complete wavelength coverage between 3 600 and 9 200 A was 
achieved, with a resolving power R = A/AA 48 000 (with 2.2 
pixels per AA resolution element), yielding a S/N of several 
hundred in V in 120, 600 and 300 sec exposures. A correspond- 
ing spectrum of /3 Ori was adopted from Commissioning II data 
(#0783, 20 sec exposure taken in November 1998). 

Data reduction was performed using the FEROS context 
in the MIDAS package (order definition, bias subtraction, sub- 
traction of scattered light, order extraction, flat-fielding, wave- 
length calibration, barycentric movement correction, merg- 
ing of the orders), as described in the FEROS documen- 
tation (http: //www.ls. eso.org/lasilla/Telescopes/2p2T/Elp5M/ 
FEROS/docu/Ferosdocu.html). Optimum extraction of the or- 
ders with cosmic ray clipping was chosen. In the spectral re- 
gion longward of ~8 900 A problems with the optimum extrac- 
tion arose due to the faintness of the signal. Standard extrac- 
tion was therefore performed in this region. No correction for 
telluric lines was made. The spectral resolution is sufficiently 
high to trace the (broadened) stellar lines in our spectrum syn- 
thesis approach even within the terrestrial O2 and H2O bands 
in many cases, see e.g. Fig. 9 of Przybilla et al. ( 2001b 1. The 
spectra were normalised by fitting a spline function to carefully 
selected continuum points. This suffices to retain the line pro- 
files of the Balmer lines in supergiants as these are rather weak 
and sampled by a single Echelle order. Finally, the spectra were 
shifted to the wavelength rest frame by accounting for a ra- 
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dial velocity v rat j as determined from cross-correlation with an 
appropriate synthetic spectrum. Excellent agreement with v ra( j 
data from the literature (see Table[2} was found. Representative 
parts of the spectra are displayed in Fig.[Q 

3. Model atmospheres 

Model atmospheres are a crucial ingredient for solving the in- 
verse problem of quantitative spectroscopy. Ideally, BA-SGs 
should be described by unified (spherically extended, hydro- 
dynamical) non-LTE atmospheres, accounting for the effects 
of line blanketing. Despite the progress made in the last three 
decades, such model atmospheres are still not available for rou- 
tine applications. In the following we investigate the suitability 
of existing models for our task, their robustness for the applica- 
tions and their limitations. The criterion applied in the end will 
be their ability to reproduce the observations in a consistent 
way. 

3.1. Comparison of contemporary model atmospheres 

Two kinds of classical (plane-parallel, hydrostatic and sta- 
tionary) model atmospheres are typically applied in the con- 
temporary literature for analyses of BA-SGs: line-blanketed 
LTE atmospheres and non-LTE H+He models (without line- 
blanketing) in radiative equilibrium. First, we discuss what ef- 
fects these different physical assumptions have on the model 
stratification and on synthetic profiles of important diagnostic 
lines. We therefore include also an LTE H+He model with- 
out line-blanketing and additionally a grey stratification in the 
comparison for two limiting cases: the least and the most lu- 
minous supergiants of our sample, 77 Leo and HD 92207, of 
luminosity class (LC) lb and Iae. The focus is on the photo- 
spheric line-formation depths, where the classical approxima- 
tions are rather appropriate - the velocities in the plasma re- 
main sub-sonic, the spatial extension of this region is small 
(only a few percent) compared to the stellar radius in most 
cases, and the BA-SGs photospheres retain their stability over 
long time scales, in contrast to their cooler progeny, the yellow 
supergiants, which are to be found in the instability strip of the 
Hertzsprung-Russell diagram, or the Luminous Blue Variables. 

The non-LTE models are computed using the code Tlusty 
(Hubeny & Lanz 1995), the LTE models are calculated with 
Atlas9 (Kurucz 1993 1, in the version of M. Lemke, as obtained 
from the CCP7 software library, and with further modifica- 
tions (Przybilla et al. 2001b) where necessary. Line blanketing 
is accounted for by using solar metallicity opacity distribution 
functions (ODFs) from Kurucz ( 1992). For the grey tempera- 
ture structure, T A — | T* s [t + qir)], exact Hopf parameters q{r) 
(Mihalas fT978l p. 72) are used. 

In Fig. |3 we compare the model atmospheres, represented 
by the temperature structure and the run of electron number 
density « e . For the less-luminous supergiant marked differences 
in the line-formation region (for metal lines, typically between 
Rosseland optical depth -1 <log tr <0) occur only between 
the line-blanketed model, which is heated due to the back- 
warming effect by <200K, and the unblanketed models. In 
particular, it appears that non-LTE effects on the atmospheric 



structure are almost negligible, reducing the local temperature 
by <50 K. The good agreement of the grey stratification and 
the unblanketed models indicates that Thomson scattering is 
largely dominating the opacity in these cases. A temperature 
rise occurs in the non-LTE model because of the recombi- 
nation of hydrogen - it is an artefact from neglecting metal 
lines. The gradient of the density rise with atmospheric depth 
is slightly flatter in the line-blanketed model. In the case of 
the highly luminous supergiant, line-blanketing effects retain 
their importance (heating by <300K), but non-LTE also be- 
comes significant (cooling by <200K). From the comparison 
with observation it is found empirically, that a model with the 
metal opacity reduced by a factor of 2 (despite the fact that 
the line analysis yields near-solar abundances) gives an over- 
all better agreement. This is slightly cooler than the model for 
solar metal abundance and may be interpreted as an empirical 
correction for unaccounted non-LTE effects on the metal blan- 
keting in the most luminous objects, i.e. ySOri and HD 92207 
in the present case. The local electron number density in the 
line-blanketed case is slightly lower than in the unblanketed 
models. We conclude from this comparison that at photospheric 
line-formation depths the importance of line blanketing out- 
weighs non-LTE effects, with the role of the latter increasing 
towards the Eddington limit, as expected. In the outermost re- 
gions the differences between the individual models become 
more pronounced. However, none of the stratifications can be 
expected to give a realistic description of the real stellar at- 
mosphere there, as spherical extension and velocity fields are 
neglected in our approach. 

The results of spectral line modelling depend on the de- 
tails of the atmospheric structure and different line strengths 
result from computations based on the different stratifications. 
Therefore, non-LTE line profiles for important stellar param- 
eter indicators are compared in Fig. [5] H6 as a representa- 
tive for the gravity-sensitive Balmer lines, a typical He 1 line, 
and features of Mg 1/11, which are commonly used for the T e g- 
determination. In the case of the Balmer and the Mg 11 lines dis- 
crimination is only possible between the line-blanketed and the 
unblanketed models. Both the LTE and non-LTE H+He model 
structure without line blanketing result in practically the same 
profiles. This changes for the lines of He 1 and Mg 1, which re- 
act sensitively to modifications of the atmospheric structure: all 
the different models lead to distinguishable line profiles. These 
differences become more pronounced at higher luminosity. 

This comparison of model structures and theoretical line 
profiles is instructive, but the choice of the best suited model 
can only be made from a confrontation with observation. It 
is shown later that the physically most sophisticated approach 
available at present, LTE with line blanketing plus non-LTE 
line formation and appropriately chosen parameters, allows 
highly consistent analyses of BA-SGs to be performed. This 
is discussed next. 

3.2. Effects of helium abundance and line blanketing 

Helium lines become visible in main sequence stars at the tran- 
sition between A- and B-types at T e g ~ 10 000 K. In supergiants 
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Fig. 2. Comparison of different model atmospheres: temperature structure (top) and electron density (bottom) for two supergiant 
models as a function of the Rosseland optical depth. On the left for the LC lb with stellar parameters according to our model for 
j] Leo, on the right for the LC Iae (HD 92207), corresponding to our sample objects at lowest and highest luminosity. Full line: 
Atlas9 line-blanketed LTE model (adopted for the analysis), dashed: Atlas9 H+He LTE model without line blanketing, dashed- 
dotted: Tlusty H+He non-LTE model without line blanketing, boxes: grey atmosphere. A second Atlas9 line-blanketed LTE 
model is displayed in the case of the Iae object (dotted), calculated with background opacities corresponding to solar metallicity. 
In our final model for HD 92207 the background opacities are reduced by an empirical factor of two in metallicity, correcting at 
least qualitatively for neglected non-LTE effects. In the inset the formation region of weak lines is displayed enlarged. 



this boundary is lowered to ~8 000 K due to stronger non-LTE 
effects in He i and by the commonly enhanced atmospheric he- 
lium abundance in these stars (see Sect.[6j. The main effect of 
an helium enhancement is the increase of the mean molecu- 
lar weight of the atmospheric material, affecting the pressure 
stratification; the decrease of the opacity and thus an effect 
on the temperature structure on the other hand is negligible 
(Kudritzki ll973> . Both effects are quantified in Fig. [4] where 
the test is performed for solar and an enhanced abundance 
of y = 0.15, while all other parameters remain fixed. The rel- 
ative increase in density strengthens with decreasing surface 
gravity. The spectrum analysis of the most luminous super- 
giants is distinctly influenced by helium enhancement, while 
farther away from the Eddington limit the effects diminish, see 
Fig. |5] Obviously, the Hei lines are notably strengthened with 
increasing abundance. For higher surface gravities, the Balmer 
lines and the ionization equilibrium of Mg i/ii are almost unaf- 



fected. On the other hand, near the Eddington limit the Balmer 
lines are noticeably broadened through an increased Stark ef- 
fect, simulating a higher surface gravity in less careful anal- 
yses. A marked strengthening of the Mg i lines is noticed, as 
the ionization balance is shifted in favour of the neutral species 
through the increased electron density. In addition, the lines 
from both ionic species of magnesium are strengthened by the 
locally increased absorber density: the combined effects result 
in a higher effective temperature from Mg i/n, when helium en- 
hancement is neglected. All studies of highly luminous BA- 
type supergiants (like a Cyg or /? Ori) to date have neglected 
atmospheric structure modifications due to helium enhance- 
ments, as these introduce an additional parameter into the anal- 
ysis. It is shown that this is not justified: a higher precision in 
the stellar parameter determination is obtained in the present 
work by explicitly accounting for this parameter. 
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Fig. 3. Comparison of theoretical non-LTE profiles for several diagnostic lines, calculated on the basis of different model at- 
mospheres for AO supergiants of LC lb (77 Leo, left) and Iae (HD 92207, right). Ordinate is normalised flux. The same line 
designations as in Fig.|2]are used, and the profiles are broadened accounting for instrumental profile, rotation and macroturbu- 
lence (see Table|2). Note that the He 1 and Mg 1 lines react sensitively to the detailed structure of the models while for the Balmer 
lines and the Mg 11 lines only the profiles of the line-blanketed and the unblanketed models can be discriminated well. 



Line blanketing is an important factor for atmospheric anal- 
yses. However, it is not only a question of whether line blanket- 
ing is considered or not, but how it is accounted for in detail. 
The two important parameters to consider here are metallici- 
ty and microturbulence, which both affect the line strengths 
and consequently the magnitude of the line blanketing effect. 
In studies of supergiants so far this has been neglected. The 
introduction of two extra parameters further complicates the 
analysis procedure and requires additional iteration steps, but 
is rewarded in terms of accuracy. 

The influence of the metallicity on the atmospheric line 
blanketing effect is displayed in Fig. 0] Here all parameters 
are kept fixed except the metallicity of the ODFs used for the 
model computations. A model sequence for four metallicities, 
spanning a range from solar to 0. 1 x solar abundances, is com- 
pared. In the LC lb supergiant model, the density structure is 
hardly affected and the local temperatures in the line-formation 
region differ by less than 100K for a change of metallicity 



within a factor of ten. This difference increases to ~200 K close 
to the Eddington limit: the higher the metallicity, the stronger 
the backwarming and the corresponding surface cooling due 
to line blanketing. In addition, the density structure is also no- 
tably altered, to a larger extent as in the case of a moderately 
increased helium abundance. Radiative acceleration diminishes 
for decreasing metallicity and thus the density rises. The corre- 
sponding effects on the line profiles are summarised in Fig. [5] 
An appreciable effect is noticed only for the highly temperature 
sensitive He 1 and Mg 1 lines at LC lb. Again, at high luminosity 
all the diagnostic lines are changed considerably, the extreme 
case being the Mg 1 line which is strengthened by a factor of 
almost three. Ignoring the metallicity effect on the line blanket- 
ing in detail will result in significantly altered stellar parame- 
ters from the analysis. 

Microturbulence has a similar impact on the line blanket- 
ing. An increase in the microturbulent velocity £ strengthens 
the backwarming effect as does an increase in the metallicity, 
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Fig. 4. Influence of the helium abundance on the atmospheric structure (bottom) and effects of metallicity (middle) and microtur- 
bulence (top set of curves) on the atmospheric line blanketing: temperature structure (upper panels) and electron density (lower 
panels). On the left the comparison is made for a LC lb object with stellar parameters matching those of our model for 77 Leo, 
on the right for the LC Iae (HD 92207). Only one parameter is modified each time: atmospheric helium abundance of y = 0.089 
(solar value, full line) and y = 0.15 (dotted line); using ODFs with [M/H] =0.0, -0.3, -0.7 and -1.0 dex (full, dotted, dashed, 
dashed-dotted lines); using ODFs with £ = 8, 4 and 2kms~ 1 (full, dotted, dashed lines). In the inset, the formation region for 
weak lines is enlarged. For convenience, the individual sets of curves have been shifted by the factors indicated. 



since a larger fraction of the radiative flux is being blocked. The 
resulting atmospheric structures from a test with ODFs at three 
different values for microturbulence, for £ = 2,4 and 8 km s , 
and otherwise unchanged parameters are shown in Fig. 0] In 
both stellar models the local temperatures in the line-formation 
region are increased by ~100K when moving from the lowest 
to the highest value of A noticeable change of the density 
structure is only seen for the LC Iae model. The correspond- 
ing changes in the line profiles are displayed in Fig. [3] Similar 
effects are found as in the case of varied metallicity, however 
they are less pronounced. 



The consistent treatment proposed here improves the sig- 
nificance of analyses, in particular for the most luminous su- 
pergiants. This also applies to the closely related line blocking, 
which is likewise treated in a consistent manner. 



3.3. Neglected effects 



Spherical extension of the stellar atmosphere is the first of a 
number of factors neglected in the current work. It becomes 
important in all cases where the atmospheric thickness is no 
longer negligible compared to the stellar radius. Observable 
quantities like the emergent flux, the colours and line equiv- 
alent widths from extended models will deviate from plane- 
parallel results for increasing extension rj (= atmospheric thick- 
ness/stellar radius), which can lead to a modified interpretation 
of the observed spectra. The expected effects on the line spec- 
trum are (mostly) reduced equivalent widths due to extra emis- 
sion from the extended outer regions and a shift in the ionisa- 
tion balance. Details on the differences between spherically ex- 
tended and plane-parallel hydrostatic LTE model atmospheres 
can be found in Fieldus et al. ( 1990 1. Note that rj will be on 
the order of a few percent for the photospheres of the objects 
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Fig. 5. Influence of the helium abundance (bottom), metallicity (middle) and microturbulence (top set of curves) on the profiles 
of diagnostic lines for our models for 77 Leo and HD 92207. The same designations as in Fig. @] are used, and the profiles are 
broadened accounting for instrumental profile, rotation and macroturbulence. Note the strong sensitivity of the lines to these 
'secondary' parameters at LC Iae. 



investigated here (adopting atmospheric thicknesses from the 
Atlas9 models and stellar radii from Table |2ji. 

Evidence for macroscopic velocity fields, i.e. a stellar wind, 
in BA-SG atmospheres is manifold, most obviously from P- 
Cygni profiles of strong lines but also from small line asym- 
metries with extra absorption in the blue wing and blue- 
shifts of the central line wavelength in less spectacular cases. 
Kudritzki ( 1992 > investigated the influence of realistic velocity 
fields from radiation driven winds on the formation of photo- 



spheric lines (in plane-parallel geometry). The subsonic out- 
flow velocity field at the base of the stellar wind strengthens 
lines that are saturated in their cores even for the moderate 
mass-loss rates typically observed for BA-SGs. Desaturation of 
the lines due to the Doppler shifts experienced by the moving 
medium is the driving mechanism for this strengthening. 

Velocity fields therefore counteract the proposed line weak- 
ening effects of sphericity, apparently leading to a close net 
cancellation for the weak line spectrum. Plane-parallel hydro- 
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static atmospheres thus seem a good approximation to spherical 
and hydrodynamic (unified) atmospheres for the modelling of 
BA-SGs, if one concentrates on the photospheric spectrum. In 
fact, first results from a comparison of classical LTE with uni- 
fied non-LTE atmospheres in the BA-SG regime indicate good 
agreement (Santolaya-Rey et al. H997t Puis et al. 120051 J. Puis, 
private communication). The situation appears to be similar in 
early B-supergiants (Dufton et al. 2005) 

BA-type supergiants have been known as photometric and 
optical spectrum variables for a long time. The most com- 
prehensive study to date in this context is that of Kaufer et 
al. J1996I IT997 >. Additional observational findings in the UV 
spectral region, in particular of the Mg 11 and Fe n resonance 
lines, are presented by Talavera & Gomez de Castro ( 1987 1 and 
Verdugo et al. dl999a> . Kaufer et al. ([1997 1 find peak-to-peak 
amplitude variations of the line strength of 29% on time scales 
of years. Unaccounted variability is therefore a potential source 
of inconsistencies when analysing data from different epochs, 
see Sect.|5]for our approach to overcome this limitation. 

BA-type supergiants are slow rotators with typical observed 
values of v sin i between 30 to 50 km s (Verdugo et al. 1999b 1. 
The modelling can therefore be treated as a 1-D problem, as 
rotationally induced oblateness of the stars is insignificant. 
Finally, magnetic fields in BA-SGs appear to be too weak to 
cause atmospheric inhomogeneities. For ySOri a weak longi- 
tudinal magnetic field, on the order of 100G was observed by 
Severny ( 1970 1. Little information on other objects is available. 

3.4. Limits of the analyses 

The spectrum synthesis technique described in the present 
work is applicable to a rather wide range of stellar parameters, 
but nevertheless it is restricted. Its scope of validity principally 
concentrates on BA-SGs and related objects of lower luminos- 
ity. This is mainly due to the limits posed by the underlying 
atmospheric models and the atomic models implemented. 

From estimates, such as that presented by Kudritzki ( 1988 
Fig. m, 9), it is inferred that non-LTE effects on the atmospheric 
structure, increasing with stellar effective temperature, will in- 
hibit analyses with the present technique of any supergiants 
above ~20000K, i.e. in the early B-types. However, main se- 
quence stars and even (sub-)giants of such spectral type are 
analysed with classical atmospheric models on a routine basis. 
For the less-luminous supergiants of mid B-type the method 
may still be applicable. This requires further investigation (in- 
cluding an extension of the model atom database to doubly- 
ionized species), but it can be expected to fail at higher lu- 
minosity. As Dufton et al. (2005 1 indicate that classical and 
unified non-LTE atmosphere analyses give similar results for 
early B-SGs, the solution for analyses of highly luminous mid- 
B supergiants will be to use non-LTE line-formation compu- 
tations based on classical line-blanketed non-LTE (instead of 
LTE) model atmospheres. 

The lower limit (in T e ff) for the applicability of the method 
is determined by several factors. At T e ff ^ 8 000 K helium lines 
disappear in the spectra of A-supergiants. Thus the helium 
abundance has to remain undetermined, introducing some un- 
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Fig. 6. Example for the occurrence of pressure inversion in 
supergiants of spectral types later than ~A4. The pressure 
stratification as a function of mass scale is displayed for two 
Atlas9 models for T e g = 7 700 K and log g = 0.42 (the com- 
puted Eddington limit, dotted line) and 0.47 (full line). The 
inset compares the respective H<5 profiles. A change in \ogg 
by 0.05 dex results in equivalent widths differing by a factor of 
2 - an artefact of the pressure bump in the line-formation depth 
and its impact on the Stark broadening. 

certainties into the analyses. Around the same temperature 
convection can be expected to set in, as shown by Simon et 
al. (2002 1 for main sequence stars, who place the boundary 
line for convection at r e ff — 8 250 K (however, no information 
is available for A-supergiants). Since the theoretical consid- 
erations of Schwarzschild dl975> . which predict only a small 
number of giant convection cells scattered over the stellar sur- 
face, atmospheric convection in supergiants has attracted little 
interest until recently. Interferometric observations of the late- 
type supergiant a Ori (Young et al. 120001 can be interpreted 
in favour of this, further strengthened by first results from 3-D 
stellar convection models (Freytag et al. l2002> . Again, no quan- 
titative information on this is available for mid and late A- 
supergiants, resulting in a potential source of systematic er- 
ror for model atmosphere and line formation computations. 
Furthermore, convective stellar envelopes give rise to chromo- 
spheres (see Dupree et al. d2005> for observational evidence 
in the cooler luminous stars), which introduce an additional 
source of UV irradiance, altering the non-LTE populations of 
the photospheric hydrogen (as compared to models without 
chromospheres) and potentially affecting the FL opacity and 
thus the stellar continuum (Przybilla & Butler 2004b 1. 

Moreover, supergiant model atmospheres cooler than T e g 
~ 8 200 K are characterised by pressure inversion (and accom- 
panying density inversion). Pressure inversion can develop 
close to the Eddington limit when radiative acceleration lo- 
cally dominates over gravity because of an opacity bump in 
the hydrogen ionization zone, which is located in the photo- 
sphere at these effective temperatures. An example of the ef- 
fect is shown in Fig. [6] where two Atlas9 models close to the 
hotter end of the pressure inversion regime are compared. A 
small change of surface gravity by ~10% results in a dras- 
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Table 1. 


Non-LTE model atoms 


Ion 


Source 


H 


Przybilla & Butler ( 2004a l 


He i 


Husfeld et al. ( 1989 1, with updated atomic data 


Ci/ii 


Przybilla et al. <2001b> 


Ni/n 


Przybilla & Butler i200fl 


Oi/a 


Przybilla et al. 1 2000 1 combined with Becker & 




Butler ( 1988 1, the latter with updated atomic data 


Mgi/n 


Przybilla et al. <2001a> 


S ii/m 


Vrancken et al. 1 1996 1, with updated atomic data 


Tin 


Becker tl998i 


Fen 


Becker 1 19981 



tic change of a factor of 2 in equivalent width of the Balmer 
lines in this (extreme) case. Despite a general strong sensitiv- 
ity of the hydrogen line equivalent widths to surface gravity 
close to the Eddington limit (see Sect. |5} this huge effect is 
an artefact of the modelling for the most part. Another pos- 
sible effect of pressure inversion on line-formation computa- 
tions concerns deviations from generally inferred trends in the 
behaviour of line strengths with stellar parameter variations. 
At slightly hotter temperatures the hydrogen line strengths de- 
crease with decreasing surface gravity. In the pressure inver- 
sion regime this trend can be compensated, and even reversed, 
by a developing pressure bump. Naturally, pressure inversion 
affects all spectral features with line-formation depths coincid- 
ing with the pressure bump because of its effect on absorber 
densities. Systematic errors for stellar parameters from ioniza- 
tion equilibria and chemical abundance determinations can be 
expected. Pressure inversion is also discussed in the context of 
hydrodynamical models (Achmad et al. 119971 Asplund 1998 1. 
It is not removed by stationary mass outflow (except for very 
high mass-loss rates not supported by observation). It is not 
initiating the stellar wind either. Abolishing the assumption of 
stationarity provides a solution to the problem as discussed by 
de Jager ( 1998 and references therein), leading to pulsations 
and enhanced mass-loss in the yellow super- and hypergiants, 
which are located in a sparsely populated region of the empiric 
Hertzsprung-Russell diagram. 

From these considerations we restrict ourselves to super- 
giant analyses at > 8 000 K using the methodology pre- 
sented here. In our opinion extension to studies of supergiants 
at cooler temperatures (of spectral types mid/late A, F, G, and 
of Cepheids) requires additional theoretical efforts in stellar 
atmosphere modelling, far beyond the scope of the present 
work. At lower luminosities the problems largely diminish so 
that the hybrid non-LTE technique provides an opportunity to 
improve on the accuracy of stellar analyses over a large and 
important part of the Hertzsprung-Russell diagram. Non-LTE 
model atoms for many of the astrophysically interesting chem- 
ical species are already available for studies of stars of later 
spectral types, as will be discussed next. 

4. Statistical equilibrium and line formation 

Detailed quantitative analyses of stellar spectra require another 
ingredient besides realistic stellar model atmospheres: an ac- 



curate modelling of the line-formation process. While our re- 
stricted modelling capacities do not allow the overall problem 
to be solved in a completely natural fashion, i.e. simultane- 
ously, experience tells us that we can split the task into sev- 
eral steps. In particular, while non-LTE effects are present in 
all cases, as photons are leaving the stellar atmosphere, they 
may be of little importance for the atmospheric structure when 
the main opacity sources remain close to LTE (H and He in 
early-type stars). Minor species - this includes trace elements 
as well as high-excitation levels of hydrogen and helium - be- 
have in this way. They can be treated in a rather coarse approx- 
imation for atmospheric structure computations, while detailed 
non-LTE calculations may be required for the analysis of their 
spectra in order to use them as stellar parameter or abundance 
indicators. We therefore chose such a hybrid approach for our 
analyses: based on line-blanketed classical model atmospheres 
we solve the statistical equilibrium and the radiative transfer 
problem for individual species in great detail. The derived level 
occupations are then used in the formal solution to calculate the 
emergent flux, considering exact line-broadening. The last two 
steps are performed using the non-LTE line-formation package 
Detail and Surface (Giddings[l98lJ Butler & Giddings 1985 1, 
which has undergone substantial extension and improvement 
over the years. In our context the inclusion of an Accelerated 
Lambda Iteration (ALI) scheme (Rybicki & Hummer 199T) is 
of primary interest, as it allows elaborate non-LTE model atoms 
to be used while keeping computational expenses moderate. 
Line blocking is taken into account via ODFs. Special care is 
required in computations for major line opacity sources, like 
iron: in order to prevent counting the line opacity twice (via 
the ODFs and as radiative transitions in the statistical equilib- 
rium computations) we adopt ODFs with a metallicity reduced 
by a factor 2. This approximately corrects for the contribution 
of the element to the total line opacity. 

At the centre of our hybrid non-LTE analysis technique 
stand sophisticated model atoms, which comprise many of 
the most important elements in the astrophysical context. An 
overview is given in Tabled summarising the references where 
the model atoms and their non-LTE behaviour were discussed 
in detail. In brief, they are characterised by the use of ac- 
curate atomic data, replacing approximate data as typically 
used in such work, by experimental data (a minority) or data 
from quantum-mechanical ab-initio computations (the bulk). 
We have profited from the efforts of the Opacity Project (OP; 
Seaton et al. 11994b and the IRON Project (IP; Hummer et 
al. Qp93 ), as well as numerous other works from the physics 
literature. In particular, a major difference to previous efforts is 
the use of large sets of accurate data for excitations via electron 
collisions. Thus, not only the (non-local) radiative processes 
driving the plasma out of LTE are treated realistically (line 
blocking is considered via Kurucz ( 1992 1 ODFs), but also the 
competing (local) processes of thermalising collisions. These 
are essential to bring line analyses from different spin systems 
of an atom/ion into agreement. A few of the older models have 
been updated/extended with respect to the original publica- 
tions. In the case of S ii/m the fits to the original photoionization 
cross-sections of the OP data (neglecting resonances due to au- 
toionizing states) have been replaced by the detailed data, and 
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Fig. 7. Comparison of non-LTE departure coefficients b, for 
Oi (top) and Fen (bottom) in HD 92207 as a function of tr. 
Displayed are the departure coefficients of the ground state 
(thick-dotted line), of the ground state of the next higher ion 
(long-dashed line) and those of the levels from which the 
diagnostic spectral lines in the optical/near-IR arise (Oi). In 
the case of Fen the departure coefficients of all non-LTE 
levels are shown. The two ionic species are representative for 
the non-LTE behaviour of the light (including the a-process 
elements) and iron group elements, see the text for further 
discussion. 



for these ions, as well as for He i and O n the features treated in 
the line-formation computations have been extended, as well as 
oscillator strengths and broadening parameters updated to more 
modern values, see Appendix 1X1 for further detail. In the fol- 
lowing we will discuss some more general conclusions, which 
can be drawn from the analysis of such a comprehensive set of 
model atoms, while referring the reader to the original publica- 
tions (see Tabled for the details of individual atoms/ions. 

There is a fundamental difference between the non-LTE be- 
haviour of the light and a-process elements on the one hand and 
the iron group elements on the other. The first group is charac- 
terised by a few valence electrons, which couple to only a few 
low-excitation states in the ground configuration that are sepa- 
rated by a large energy gap (several eV) from the higher excited 
levels. These excited states in turn show a similar structure on 
a smaller scale: a few (pseudo-)metastable levels are detached 



from the remainder by a ~2eV gap, see e.g. the Grotrian dia- 
grams in the references given in Tabled Collisional processes 
are effective in coupling levels either below or above the en- 
ergy gaps, such that these are in or close to LTE relative to 
each other. The levels of highest excitation couple to the ground 
state of the next higher ionization stage via collisions. On the 
other hand, only a few electrons in the high-velocity tail of the 
Maxwell distribution are energetic enough in the atmospheres 
of BA-type stars to pass the first gap, and considerably more, 
but a minority nonetheless, the second gap. This favours strong 
non-LTE overpopulations of the excited (pseudo-)metastable 
states, as shown exemplarily for Oi in Fig. manifested in 
non-LTE departure coefficients bj = nf urE /nf rE > 1, where the 
rii are non-LTE and LTE level occupation numbers, respec- 
tively. Therefore, the diagnostic lines in the optical/near-IR 
from the light and a-process elements are typically subject to 
non-LTE strengthening. 

On the other hand, the electrons in the open 3c/-shell of the 
iron group elements give rise to numerous energetically close 
levels throughout the whole atomic structure. Photoionizations 
from the ground states of the single-ionized iron group ele- 
ments are typically not very effective, as the ionization thresh- 
olds fall short of the Lyman limit where the stellar flux is neg- 
ligible. The situation is different for the well populated lev- 
els a few eV above the ground state. They show the largest 
non-LTE depopulations (see Fig. for the example of Fen). 
Because of the collisional coupling the ground states also be- 
come depopulated, but to a lower degree. Nonetheless, a net 
overionization is established, resulting in a non-LTE overpop- 
ulation of the main ionization stage (Fein in Fig. 0- note 
that no Fe in lines are observed in the optical/near-IR for this 
star). Again, the highest excitation levels of the lower ioniza- 
tion stage couple collisionally to this and are thus also over- 
populated. A continuous distribution of departure coefficients 
results, with the b, of the upper levels of the transitions typi- 
cally larger than those of the lower levels. Consequently, the 
optical/near-IR lines from the iron group elements experience 
non-LTE weakening. Therefore, LTE analyses of luminous BA- 
SGs will tend to systematically overestimate abundances of the 
light and a-process elements, and to underestimate abundances 
of the iron group elements. These effects will be quantified in 
Sect. |6]where the comparison of non-LTE and LTE computa- 
tions with observation is made. Note however that non-LTE 
computations, when performed with inaccurate atomic data, 
also bear the risk of introducing systematic errors to the analy- 
sis, see e.g. the discussions in Przybilla & Butler (2001 2004a I. 
This is because of the nature of statistical equilibrium, where 
(de-)population mechanisms couple all energy levels with each 
other. The predictive power of non-LTE computations is there- 
fore only as good as the models atoms used. 

The spectral lines treated in our non-LTE approach com- 
prise about -70% of the observed optical/near-IR features in 
BA-SGs. This includes in particular most of those of large and 
intermediate line strength. The remainder of the observed lines 
is typically weak, except for several Sin and Cm transitions. 
In order to achieve (near) completeness we incorporate these 
and another ten chemical species into our spectrum synthesis in 
LTE. It can be expected that this approach will introduce some 



12 



N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



systematic deficiencies, in particular for the most luminous ob- 
jects, as will be indicated by the analysis in Sect. [6] However, 
the solution for high-resolution observations is to draw no fur- 
ther conclusions from these species. Interpretation of medium- 
resolution spectra (see Sect. [8ji on the other hand would suf- 
fer more from their absence than from their less accurate treat- 
ment, as they typically contribute only small blends to the main 
diagnostic features. 

The currently used line lists comprise several ten-thousand 
transitions, covering the classical blue region for spectral anal- 
yses between ~4000 and 5 000 A well. At longer wavelengths 
a number of features are missing in our spectrum synthesis 
computations, mostly lines from highly-excited levels of the 
iron group elements. However, these are intrinsically weak 
(with equivalent widths Wa < 10 mA) and typically isolated, so 
that their absence will hardly be noticed when compared to 
intermediate-resolution observations. 

Where it was once a supercomputing application, compre- 
hensive non-LTE modelling like the present can now be made 
on workstations or PCs. Typical running times to achieve con- 
vergence in the statistical equilibrium and radiative transfer cal- 
culations with Detail range from ~10min for the most simple 
model atoms to a couple of hours for the Fe n model on a 3 GHz 
P4 CPU - for one set of parameters. The formal solution for a 
total of 2-3 xlO 5 frequency points with Surface requires ~2Q- 
30 CPU min. Because of the highly iterative and interactive na- 
ture of our analysis procedure (see next two sections) the total 
time for a comprehensive study of one high-resolution, high- 
S/N spectrum, assuming excellent wavelength coverage from 
around the Balmer jump to ~9 000A as typically achieved by 
modern Echelle spectrographs, amounts to 1-2 weeks for the 
experienced user. This is the price to pay for overcoming the 
restrictions of contemporary BA-SG abundance studies and ad- 
vancing them from factor 2-3 accuracy astronomy to precision 
astrophysics. 

Future extensions of the present work will aim to provide 
more non-LTE model atoms, but these quite often require the 
necessary atomic data to be calculated first, as many data are 
still unavailable in present-day literature. In particular, the sta- 
tus of collisional data has largely to be improved. Note that the 
Si ii ion in the silicon model atom by Becker & Butler ( 1990) 
is only rudimentary (though sufficient for their main purpose, 
analyses of early B-stars). Several energy levels involved in the 
observed transitions in BA-SGs are missing, and test calcula- 
tions indicate a wide spread of the non-LTE abundances from 
the remaining transitions. On the basis of these findings we re- 
frain from using the model atom for analyses in the present 
work, but wish to emphasise that qualitatively the correct non- 
LTE behaviour - line strengthening - is predicted. 

5. Determination of stellar parameters 

Stellar parameter determination for BA-SGs is complicated be- 
cause of the intrinsic variability of these objects. Observational 
data from different epochs may therefore reflect different phys- 
ical states of the stellar atmosphere, unless these are obtained 
(nearly) simultaneously. While we can expect the mass of a 
supergiant and its surface abundances to be conserved on hu- 



man timescales, other stellar parameters may not. The observed 
variability patterns of light curves, radial velocities and line 
profiles can be interpreted in favour of a mix of radial and non- 
radial pulsations (see e.g. Kaufer et al. 1 19971 1. which may af- 
fect T e g and log g, that in turn determine colours, the spectral 
energy distribution and the spectra. Even the stellar luminos- 
ity may be subject to small changes (Dorfi & Gautschy 2000 1. 
The parameter variations are not as pronounced as in other 
variable stars, but use of information from different epochs 
can potentially introduce systematic uncertainties into high- 
precision analyses. A solution is to derive the desired quantities 
from one set of observational data. Echelle spectra are the best 
choice, containing all information required for the derivation of 
stellar parameters and elemental abundances. With few excep- 
tions, use of photometric data (from other sources) can thus be 
avoided. In the following we discuss the details of our spectro- 
scopic approach for the parameter determination of BA-SGs. 
Spectrophotometry is only briefly considered for consistency 
checks. Finally, our findings are compared with previous anal- 
yses of the two well studied standards 77 Leo and/30ri. 

5.1. Spectroscopic indicators 

The hydrogen lines are the most noticeable features in the spec- 
tra of BA-type stars. Their broadening by the linear Stark effect 
gives a sensitive surface gravity indicator in the mostly ionized 
atmospheric plasma. We employ the recent Stark broadening 
tables of Stehle & Hutcheon (1999 SH), which compare well 
with the classically used data from Vidal et al. J1973I VCS; 
with extensions of the grids by Schoning & Butler, private com- 
munication) in the BA-SG regime. The SH tables have not only 
the advantage of being based on the more sophisticated theory 
but they also cover the Paschen (and Lyman) series for transi- 
tions up to principal quantum number n = 30. A detailed dis- 
cussion of our hydrogen non-LTE line-formation computations 
for BA-SGs can be found in Przybilla & Butler (2004a 1, where 
lines of the Brackett and Pfund series were also studied. In 
short, it has been shown that except for the lower series mem- 
bers - which are affected or even dominated by the stellar wind 
- excellent consistency can be achieved from all available indi- 
cators, if accurate data for excitation via electron collisions are 
accounted for. Note that the only other study of H lines from 
these four series, in the prototype A-SG Deneb (Aufdenberg et 
al. |2002i, fails in this. Their non-LTE model atmosphere and 
line-formation computations produce much too strong near-IR 
features because of inaccurate collisional data in their hydro- 
gen model atom, as indicated by our findings. The problems 
are resolved using the improved effective collision strengths 
(P. Hauschildt, private communication). 

Our method for deriving surface gravities from the hydro- 
gen lines deviates in a few details from the usual approach, 
which is based on the modelling of the Hy and H6 line wings. 
The effect of variations of log g and T e g on the H6 profiles of 
77 Leo and HD 92207 are displayed in Fig. [8] We have chosen 
parameter offsets to our finally adopted values which we con- 
sider conservative estimates for the uncertainties of our anal- 
ysis. It is clear that the internal accuracy of our modelling is 
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Fig. 8. Impact of stellar parameter variations on non-LTE line profile fits for H6 in rj Leo (left) and HD 92207 (right). Spectrum 
synthesis for the adopted parameters (see Table|2] thick red line) and varied parameters, as indicated (thin lines), are compared to 
observation. A vertical shift by 0.5 units has been applied to the upper profiles. Note the decrease of the Balmer line strength in 
the progression from LC lb (77 Leo) to Iae (HD 92207), and the increased sensitivity to the surface gravity parameter. In fact, the 
comparison indicates that for the more luminous supergiants uncertainties in the surface gravity determination become as low as 
0.05 dex (dotted lines). The sensitivity to reasonable r e ff-changes is almost negligible. Wind emission slightly contaminates the 
red wing of H6 in HD 92207. 
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Fig. 9. Temperature determination for 77 Leo using the Mg i/n non-LTE ionization equilibrium. Displayed are the observed line 
profiles for some of the strategic lines, Mgi A 5 183 (left panel) and Mgii /14390 and Hei /14387 (right panel), and the best 
non-LTE fit for stellar parameters as given in Table|2(thick red line); theoretical profiles for varied parameters are also shown 
(thin lines, as indicated). A vertical shift by 0.3 units has been applied to the upper set of profiles. Note the strong sensitivity of 
the minor ionic species, Mg 1, to parameter changes, while Mg 11 is virtually unaffected. 



higher, amounting to < 0.10 dex at LC lb and as accurate as 
0.05 dex at LC Iae. The temperature sensitivity is smaller, but 
non-negligible. Note the strong decrease of the H6 line strength 
in the luminosity progression. Close to the Eddington limit a 
modification by 0.05 dex in logg can amount to a change in 
W,\ of the Balmer lines by 15%. Note also that in the less lu- 
minous supergiant not only the line wings but also the entire 
line profile and thus also the equivalent width is accurately 
reproduced (which becomes important for intermediate- and 
low-resolution studies), indicating that the model atmosphere 
is sufficiently realistic over all line-formation depths. At the 
highest luminosities the accuracy of the spectrum synthesis de- 
grades somewhat, because of stellar wind and sphericity ef- 



fects. However, excellent agreement can be restored even for 
the higher Balmer lines, and for the accessible Paschen series 
members (see Fig. 11 in Przybilla & Butler 2004a 1, which are 
formed even deeper in the atmosphere. This also includes good 
reproduction of the series limits and their transition into the 
continua. Moreover, we also account for the effects of non-solar 
helium abundances on the density structure (see Sect. 13.21 . 
which is mandatory for achieving consistent results, but which 
has been ignored in all BA-SG analyses so far. 

In order to resolve the ambiguity in r e ff/logg another in- 
dicator has to be applied. In a spectroscopic approach this is 
typically the ionization equilibrium of one chemical species. 
Potential ionization equilibria useful for r e ff and log g determi- 
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9200 9400 9600 9800 10000 1 1500 12000 12500 

Fig. 10. Fit diagrams of temperature- and gravity-sensitive indicators for rj Leo and /3 Ori in the r e ff-log g-plane. The curves are 
parameterised by surface helium abundance y. The intersection determines T e s and logg, with uncertainties as indicated. 



nations in optical/near-IR spectroscopy of BA-SGs are C i/ii, 
Ni/ii, Oi/ii, Mgi/n, Ali/n, Aln/m, Sin/m, Sn/m, Fei/n and 
Fe ii/iii. Of these we can use only half because model atoms 
are currently unavailable for the rest. The lines of the minor 
ionic species are highly sensitive to temperature and electron 
density changes, while the weak lines of the major ionic species 
are excellent abundance indicators. A suitable set of parameters 
T e sl log g is found in the case where both ionic species indicate 
the same elemental abundance within the individual error mar- 
gins. However, great care has to be practised in the modelling, 
using carefully selected model atmospheres (see Sect. [3 and 
sophisticated non-LTE techniques (see Sects. 0and|6jl. 

An example is displayed in Fig. [9] where tests on the Mg i/ii 
ionization balance in one of the sample objects are performed. 
Results from the best fit obtained in the detailed non-LTE 
analysis are compared with those from conservative parame- 
ter studies at an unchanged elemental abundance. The sensitiv- 
ity of the predicted Mg i line strengths to changes of T e ff and 
log g within the error margins is high, indicating internal un- 
certainties of Ar e ff < 100 K and of A log g < 0. 10 dex. The pro- 
file changes by conservative parameter variation are similar to 
those achieved by abundance variations of the order 0. 1 dex. 
On the other hand, no perceptible changes are seen in the Mg n 
line for the same parameter variations. All Mgi/n lines behave 
in a similar way, and the same qualitative characteristics are 
shared with ionization equilibria from other elements. 

Both, hydrogen profiles and ionization equilibria are af- 
fected by the helium abundance, which therefore has to be si- 
multaneously constrained as an additional parameter. Enhanced 
helium abundances manifest in broadened hydrogen lines and 
shifts in ionization equilibria - in BA-SGs the spectral lines of 
the minor ionization species are strengthened. Thus, parameter 
studies assuming a solar helium abundance will tend to derive 
systematically higher effective temperatures and surface grav- 
ities. In the BA-SGs the transitions of Hei are typically weak. 
The helium abundance y (by number) can therefore be inferred 
from line profile fits to all available features. Typical uncertain- 
ties in the determination of stellar helium abundances amount 
to ~10% because of the availability of excellent atomic data. 



The whole process of the basic atmospheric parameter de- 
termination can be summarised in an r e ff-log g diagram, which 
is done exemplarily for rj Leo and /30n in Fig.E3 The inter- 
section of the different loci marks the appropriate set of stel- 
lar parameters. Note that these diagnostic diagrams differ from 
those typically found in the literature in two respects: first, the 
different indicators lead to consistent stellar parameters, and 
second, the loci are parameterised with the helium abundance, 
a necessary procedure to obtain the excellent agreement. When 
deriving the stellar parameters one can profit from the slow and 
rather predictable behaviour of the derived helium abundance 
with varying r e( f/log g. 

The microturbulent velocity £ is determined in the usual 
way by forcing elemental abundances to be independent of 
equivalent widths. The difference to previous studies of BA- 
SGs is that this is done on the basis of non-LTE line-formation 
(see Fig. 1131 and the discussion in Sect. [6). An LTE analysis 
would tend to find higher values in order to compensate non- 
LTE line strengthening. Typically, the extensive iron spectrum 
is used for the microturbulence analysis. Because the Fe n non- 
LTE computations are by far the most time-intensive, alterna- 
tives have to be considered here. The Tin and in particular the 
N i spectra are excellent replacements, the latter distinguished 
by more accurate atomic data. Later, the iron lines can be used 
to verify the ^-determination. Note that this is the only occa- 
sion where equivalent widths are used in the analysis process - 
elsewhere line profile fits are preferred. 

Good starting estimates for the microturbulent velocity in 
BA-SGs are £ = 4, 6 and 8kms 1 for objects of LC lb, lab 
and la, respectively. After the determination of an improved 
value of £ the model atmosphere has to be recalculated in some 
cases, and small corrections to r e jf/log g/y may become appli- 
cable. Only one iteration step is typically necessary to reestab- 
lish consistency in the atmospheric parameters. The uncertain- 
ties amount to typically ±lkms _1 . In the present study the 
different microturbulence indicators give a single value for £ 
within these error margins. This is consistent with findings 
from recent (LTE) work on less-luminous A-type supergiants 
(Venn ll995al [T999l. while older studies based on more simple 
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Table 2. Basic properties and stellar parameters of the sample 
stars 





HD 87737 


HD 111613 


HD 92207 


HD 34085 


Name 


^Leo 






B Ori, Rigel 


Association 11 


Field 


Cen OBI 


Car OBI 


Ori OB 1 


Spectral Type b 


AO lb 


A2 Iabe 


AOIae 


B8 Iae: 


a (J2000) b 


1007 19.95 


1251 17.98 


10 37 27.07 


05 14 32.27 


6 (J2000) b 


+ 1645 45.6 


-60 19 47.2 


-58 44 00.0 


-08 1205.9 


1 K > 


?1Q 






209.24 


b(°) b 


+50.75 


+2.54 


-0.26 


-25.25 


n (mas) c 


1.53+0.77 


1.09+0.62 


0.40+0.53 


4.22+0.81 


rf(pc) d 


630±90 


2 290+220 


3 020+290 


360+40 


R g (kpc) 


8.25+0.83 


6.97+0.84 


7.66+0.77 


8.23+0.82 


9 D (mas) e 








2.55+0.05 


v rad (kms _1 ) b 


+3.3+0.9 


-21.0+2.0 


-8.5+5.0 


+20.7+0.9 


jx a (masyr') c 


-1.94+0.92 


-5.26+0.55 


-7.46+0.53 


1.87+0.77 


He (mas yr') c 


-0.53+0.43 


-1.09+0.39 


3.11+0.44 


-0.56+0.49 


U(kms- [ ) 


4 


-49 


-103 


-6 


V(km s~') 


222 


211 


201 


214 


W(kms~') 


6 


-6 


-6 


1 


Atmospheric: 










r e ff(K) 


9600+150 


9150+150 


9 500+200 


12000+200 


logg(cgs) 


2.00+0.15 


1.45+0.10 


1.20+0.10 


1.75+0.10 


y 


0.13+0.02 0.105+0.02 


0.12+0.02 


0.135+0.02 


[M/H] (dex) 


-0.04+0.03 


-0.11+0.03 


-0.09+0.07 


-0.06+0.10 


f(kms _l ) 


4+1 


7+1 


8+1 


7+1 


((kms" 1 ) 


16+2 


21+3 


20+5 


22+5 


v sin i (km s -1 ) 


0+3 


19+3 


30+5 


36+5 


Photometric: 










V(mag) f 


3.52 


5.72 


5.45 


0.12 


B-V f 


-0.03 


+0.38 


+0.50 


-0.03 


U -B f 


-0.21 


-0.10 


-0.24 


-0.66 


E(B - V) 


0.02 


0.39 


0.48 8 


0.05 


(m - M) d 


9.0+0.3 


11.8+0.2 


12.4+0.2 


7.8+0.2 


My 


-5.54+0.3 


-7.29+0.2 


-8.82+0.2 


-7.84+0.2 


B.C. 


-0.29 


-0.23 


-0.34 


-0.78 


M b0] 


-5.83+0.3 


-7.52+0.2 


-9.16+0.2 


-8.62+0.2 



Physical: 
logL/L 

M/M ZAMS 

M/M| vo1 
M/M^ cc 
Tevoi (Myr) 



4.23+0.12 4.90+0.08 5.56+0.08 5.34+0.08 



47+8 
10+1 
10+1 
8+4 
25+5 



112+12 
16+1 
15+1 
13+4 
14+2 



223+24 
30+3 
25+3 
29+10 
7+1 



109+12 
24+3 
21+3 
24+8 
8+1 



a Blaha & Humphreys i 19891 b adopted from the Simbad database 
at CDS c Perryman et al. 119971 d see text c Hanbury Brown et 
al. f Nicolet IT97BI 8 with R v = 3.9. see Sect.l5~3l 

model atmospheres and less accurate oscillator strengths had 
to invoke different values for various elemental species, or a 
depth-dependent £ (e.g. Rosendhal ll970l Avdin [T972l . 

Microturbulence also has to be considered as additional 
non-thermal broadening agent in the radiative transfer and sta- 
tistical equilibrium computations. The Doppler width is then 
given by A/Id = 'Wc(v? + £ 2 )'^ 2 , where Aq is the rest wave- 
length of the transition, c the speed of light and v t h the thermal 
velocity of the chemical species of interest. The main effect is 
a broadening of the frequency bandwidth for absorption in as- 
sociation with a shift in line-formation depth, typically leading 



to a net strengthening of individual lines by different amounts, 
see e.g. Przybilla et al. J2000II200 1 all200 1 bl for further details. 

Individual metal abundances are typically of secondary im- 
portance for the computation of model atmospheres, as elemen- 
tal ratios are remarkably constant (i.e. close to solar) over a 
wide variety of stars, except for the chemically peculiar. We 
determine the stellar metallicity [M/H] 1 , which is the more im- 
portant parameter, as the arithmetic mean from five elements, 
for which non-LTE computations can be done and which are 
unaffected by mixing processes: [M/H] = ([O/H] + [Mg/H] + 
[S/H] + [Ti/H] + [Fe/H])/5. High weight is thus given to the a- 
process elements, which have a different nucleosynthesis his- 
tory than the iron group elements, but [a/Fe] ~ can be ex- 
pected for these Population I objects. Because of the impor- 
tance of metallicity for line blanketing a further iteration step 
may be required in the parameter determination. 

Finally, the projected rotational velocity vsin/ and the 
(radial-tangential) macroturbulent velocity ( are derived from 
a comparison of observed line profiles with the spectrum syn- 
thesis. Single transitions as well as line blends should be 
used for this, as they contain some complementary informa- 
tion. Both quantities are treated as free parameters to obtain 
a best fit via convolution of the synthetic profile with rota- 
tion and macroturbulence profiles (Gray 1992ai, while also ac- 
counting for a Gaussian instrumental profile. An example is 
shown in Fig. \^\ where the comparison with observation in- 
dicates pure macroturbulence broadening for rj Leo. The the- 
oretical profiles intersect the observed profiles if the spectral 
lines are broadened by rotation alone, resulting in slightly too 
broad line cores and insufficiently broad line wings. Typically, 
both parameters are non-zero, with the macroturbulent velocity 
amounting to less than twice the sonic velocity in the atmo- 
spheric plasma. Macroturbulence is suggested to be related to 
surface motions caused by (high-order) nonradial oscillations 
(Lucy 1976 1. Weak lines should be used for the vsin/ and £- 
determinations in supergiants in order to avoid systematic un- 
certainties due to asymmetries introduced by the outflowing ve- 
locity field in the strong lines. 

5.2. Stellar parameters of the sample objects 

We summarise the basic properties and derived stellar parame- 
ters of the sample supergiants in Table|2] The first block of in- 
formation concentrates on quantities which in principle can be 
deduced directly, with only little modelling involved. Besides 
Henry-Draper catalogue numbers, alternative names, informa- 
tion on association membership and spectral classification and 
equatorial and galactic coordinates are given. Hipparcos par- 
allaxes it have been included for completeness, as none of the 
measurements is of sufficient statistical significance. The stel- 
lar distances d are therefore deduced by other means. 

A distance modulus of 8™ 5 (viz. 500 pc) is found for the 



Ori OB 1 association by Blaha & Humphreys ( 1989). However, 
/3 Ori shows a smaller radial velocity than the mean of the Ori 



1 using the usual logarithmic notations [X] = log e(X) il - log e(X) q , 
with \ogE(X) = log(N x /N H ) + 12, the e{X) being the abundance of 
element X and the N, number densities 
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OBI association (~30kms _1 ). Therefore, the association dis- 
tance gives only an upper limit, if one assumes that the star 
was formed near the centre of Ori OBI. During its lifetime 
jSOri could have moved ~100pc from its formation site. The 
finally adopted distance of 360 pc is indicated by Hoffieit & 
Jaschek ( 1982 ), who associate B Ori with the t Ori Rl complex. 

Recent studies of the cluster NGC4755 found a distance 
modulus of 11.6+0.2 mag, with a mean value of E(B-V) of 
0.41+0.05 mag (Sagar & Cannon [1995} and 0.36 ±0.02 mag 
(Sanneret al. 2001 1. However, HD 1 1 1613 is situated at a rather 
large distance from the cluster centre and was not observed in 
either study. In a previous work (Dachs & Kaiser 1984} the 
object was found to be slightly behind the cluster by 0.2 mag. 
Consequently, this difference is accounted for in the present 
study, otherwise using the modern distance. 

The line of sight towards the Car OBI association coin- 
cides with a Galactic spiral arm, such that the star population 
in Car OB 1 is distributed in depth over 2 to 3 kpc (Shobbrook 
& Lynga 1994 1. A more decisive constraint on the distance of 
HD 92207 is indicated by Carraro et al. (2001 ), who argue that 
the star may be associated with the cluster NGC 3324. The star 
will not be gravitationally bound because of its large proper 
motion, but appears to be spatially close to the cluster, and may 
have been ejected. We therefore adopt the cluster distance mod- 
ulus for HD 92207 with an increased error margin. Finally, the 
distance to the field star rj Leo can only be estimated, based on 
its spectroscopic parallax. 

Galactocentric distances of the stars are calculated from 
the coordinates and d, using a galactocentric solar distance of 
Ro = 7.94+0.42 kpc (Eisenhauer et al. 2003 1. In one case an 
interferometrically determined true angular diameter 0u (al- 
lowing for limb darkening) is available from the literature. 
Our measured radial velocities v ra( j (from cross-correlation with 
synthetic spectra) are compatible with the literature values, 
which in combination with proper motions fi a and fi$ are 
used to calculate galactocentric velocities U, V and W, as- 
suming standard values for the solar motion (U = 1 0.00 km s~ 1 , 
V = 5.23kms- 1 ,W = 7.17kms- 1 , Dehnen & Binney 199$ rel- 
ative to the local standard of rest (220 km s -1 , Kerr & Lynden- 
Bell lT986Tl . 

In the second block our results from the spectroscopic stel- 
lar parameter determination are summarised. Photometric data 
are collected in the third block. Observed visual magnitudes 
and colours in the Johnson system are used to derive the colour 
excess E(B - V) by comparison with synthetic colours. Note 
that for HD 111613 and HD 92207 the derived colour excess 
is in excellent agreement with literature values for their par- 
ent clusters. From the distance moduli (m - M)o absolute vi- 
sual magnitudes My are calculated, which are transferred to 
absolute bolometric magnitudes Mboi by application of bolo- 
metric corrections B.C. from the model atmosphere compu- 
tations. Note also our previous comments on the photometric 
variability of BA-SGs, which can introduce some systematic 
uncertainty into our analysis, as the photometry and our spectra 
may reflect slightly different physical states of the stellar atmo- 
spheres. The sensitivity of the Johnson colours to r e ff-changes 
is not high enough to use them as an alternative indicator for a 
precise determination of the stellar effective temperature. 



Finally, we derive the physical parameters luminosity L, 
stellar radius R and spectroscopic mass M spec of the sam- 
ple stars by combining atmospheric parameters and photom- 
etry. Thus the determined stellar radius of BOn is in good 
agreement with that derived from the angular diameter mea- 
surement (99+llR Q ). From comparison with stellar evolu- 
tion computations (Meynet & Maeder l2003t zero-age main- 
sequence masses M ZAMS and evolutionary masses M evo1 are 
determined, and the evolutionary age r eV oi- The stars were of 
early-B and late-O spectral type on the main sequence. The re- 
sults will be discussed in detail in Sect. Note that the rather 
large uncertainties in the distance determination dominate the 
error budget of the physical parameters. 

5.3. Spectrophotometry 

The principal aim of the present study is a self-consistent 
method for the analysis of BA-type supergiants in the whole. 
This requires above all the reproduction of the stellar spec- 
tral energy distribution (SED). A comparison of our model 
fluxes for three of the sample supergiants with observation is 
made in Fig. [2] The observational database consists of IUE 
spectrophotometry (as obtained from the IUE Final Archive) 
and photometry in several passbands (broad-band to small- 
band) from the near-UV to near-IR, in the Johnson (Morel 
& Magn enat IT9781 Ducati 120021 for HD 9220 7), Ge neva 
(Rufener 1988 1, Stromgren (Hauck & Mermilliod |T998 ) and 
Walraven systems (de Geus et al. 1990). We have omitted 
HD 111613 from the comparison because crucial UV spec- 
trophotometry is unavailable for this star. The observations 
have been de-reddened using a reddening law according to 
Cardelli et al. ( 1989 1. Overall, good to excellent agreement is 
found. This verifies our spectroscopic stellar parameter deter- 
mination. 

Use of SED information can provide an additional T e ff- 
indicator, independent of weak spectral lines used in the ion- 
ization equilibria approach. This makes SED-fitting attractive 
for extragalactic applications, as it has the potential to improve 
on the modelling situation when only intermediate-resolution 
spectra are available (see Sect. [8}. As the present work concen- 
trates on the spectroscopic approach to BA-SGs analyses, we 
will report on a detailed investigation of this elsewhere. 

5.4. Comparison with previous analyses 

The two MK standards 77 Leo and B Ori have been analysed 
with model atmosphere techniques before. A comparison of 
our derived stellar parameters with literature values is made in 
Table|3] The methods used for the parameter determination are 
indicated and comments on the model atmospheres and the ob- 
servational data are given. For the comparison we have omitted 
publications earlier than the 1970ies. 

Good agreement of the present parameters for rj Leo with 
those of Venn ( 1995ai is found, as the methods employed are 
comparable. Other authors find a substantially hotter T & g for 
this star, while the values for surface gravity and microturbu- 
lence are comparable. All these values are based on less elab- 
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Table 3. Comparison of stellar parameters of rj Leo and fl Ori from the literature 



Source 


r eff (K) 


log^(cgs) 


— — . T 

f (kms -1 


) Method 


Notes 










HF) 87737 




This work 


n ^nn_i_ 1 zr\ 
y oUU± 1 jU 


Z.UU+U.l J 


4+1 


"\TT TC IT T VT t/it \ /I , /,, ff~* I / TT \ 

JNL.1C, Ml, JN I/II, Jvlg I/II, (L, I/II), 












spectrophotometry 






V /UU+ZUU 


7 n_i_n 7 


A _i_ 1 
4± 1 


riy, lNL,1 H JVlg I/II 


Kurucz < 1 993 ► atmospheres 


LOuei ei ai. i lyyz \ 


1U ZUU± J IV 


1 Qj_n a 
1 .y+u.^ 


^ A-ufl 7 
J.4±U. / 


T TP P^ t/tt 
J_, rc I/II 


Kurucz < 1979 ► atmospheres, 












equivalent widths from Wolf 119711 


L,dmDen ei ai. 1 1700 1 


1 n son 
iu jyjv 


z.z 




otromgren photometry + rip 


Kurucz \ 1 979 1 atmospheres 


Wr»lf J 1 Q7 1 1 
WOII / 1| 


1U JUU 


Z.UO±U.ZU 


z. . . iu 


Balmer lines, Balmer jump, 


early LTE atmospheres, no line- 










T TF Mp- t/tt Fp t/tt 

LJ Li IVlg 1/11, 1 fc> J/ 11 


hlantptintr ntintrKyrMntiip cnpptrn 










HD 34085 




This work 


12000±200 


1.75+0.10 


7±1 


NLTE H 1, N 1/11, 1/11, S ii/m, 












spectrophotometry 




McErlean et al. fTW)\ 


13 000±1000 


1.75+0.2 


10 


NLTE Hy, H5, Sin/m 


Hubeny fPMl H+He NLTE models, 












no line-blanketing 


Israelian et al. 119971 


13 000±500 


1.6±0.1 


7 


NLTE Hy, H<5, Si 


Hubeny fTWgl H+He NLTE models, 












no line-blanketing 


Takeda 119941 


13 000±500 


2.0±0.3 


7 


LTE Hy, H<5, spectrophotometry 


Kurucz 1 1979 1 atmospheres 


Stalio et al. fPJTl 


12070±160 






angular diameter, total flux 


early line-blanketed LTE atmospheres 



orate model atmospheres, with less or no line-blanketing, and 
they completely ignore non-LTE effects on the line formation. 
Note that a T e g of 9380 K is obtained from Stromgren photom- 
etry, using the calibration of Gray ( 1992b I. 

For jSOri basically two disjunct values for the effec- 
tive temperature are found in the literature. Model atmo- 
sphere analyses so far found a systematically higher T e g, 
which is likely to be due in good part to ignored/reduced 
line-blanketing. The measurements of surface gravity and 
microturbulent velocity are again similar. The other group 
of T e g analyses mostly find significantly lower values than 
that derived in the present work, from measured total fluxes 
and interferometric radius determinations, or from the in- 
frared flux m ethod: 1 1 800 ± 300/1 1 70 ± 20 (Nandy & 
Schmidt 119751, 11 550+170 (Code et al. 119751, 11 4 10 ±330 
(Beeckmans IT9771, 1 1780 (Underhill et al. 119791, 11 014 
(Blackwell et al. fT980l . 1 1 3 80 (Un derhill & Doazan fT982l and 
1 1023/11 453 K (Glushneva [19851 . These methods are prone 
to systematic errors from inappropriate corrections for inter- 
stellar absorption. The authors all supposed zero interstellar 
extinction. In fact, the only such study that accounts for a 
non-zero E(B - V) (Stalio et al. H977l +0.04 vs. +0.05 as de- 
rived here) finds a temperature in excellent agreement with the 
present value, see Table [3] However, using LTE Fe ii/m and 
Sin/in ionization equilibria and Balmer profiles they derive 
r eff >13 000K. 



N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



19 



Fig. 12.Computed 

(red 
full 
lines) 
and 
mea- 
sured 
spec- 
tral 
en- 
ergy 
dis- 
tri- 
bu- 
tions 
for 
the 
sam- 
ple 
stars 
from 
the 
UV 
to 
the 
near- 
IR 
in 
the 
K- 

band. 
Overall, 
good 
to 
ex- 
cel- 
lent 
agree- 
ment 
is 

ob- 
tained. 
Displayed 
are 
IUE 
spec- 
tropho- 
tom- 
e- 
try 
(full 
lines), 
Johnson 
(boxes), 
Stromgren 
(tri- 
an- 
gles), 
Walraven 
(filled 
cir- 
cles) 
and 

Geneva 



Table 4. Elemental abundances in the sample stars 



Element 


Sun a 


jjheo 


T T~T"\ 1 1 i /II 

HD 111613 


HD 92207 


/JOri 


Gal AI b 


Gal BV C 


Gal BV d 


Gal BV C 


Gal BV* 


He i 


10.99±0.02 


11.18+0.04 (14) 


7i. 07+0.05 (10) 


11.14+0. 04 (10) 


11.19+0.04 (15) 






11.13±0.22 


11.05±0.10 




C i 


8.52±0.06 


7.94±0. 10(A) 


8.10+0.07 (J) 






Ci 7 A . S\ 

8.14+0.13 










C ii 


8.52±0.06 


8.10±0.09 (3) 


8.24+0.06 (3) 


8.33 (1) 


5.i5±0.05(3) 




8.27±0.13 


8.20±0.10 


8.22±0.15 


7.87+0.16 


N i 


7.92±0.06 


8.41±0.09 (20) 


8.40+0.10 '(16) 


5.25+0.04(11) 


5.50+0.07(11) 


8.37±0.21* 










Nil 


7.92±0.06 


8.32 (1) 


<5.i6(l) 


8.25(1) 


8.51+0.06 (16) 




7.63±0.15 


7.75±0.27 


7.78±0.27 


7.90+0.22 


Oi 


8.83±0.06 


8. 78+0.05 (13) 


& 70+0.04 (9) 


8.79±0.07 (6) 


8.78+0.04 (is) 


8.77+0.12 










O ii 


8.83±0.06 








8.83±0.03 (5) 




8.55±0.14 


8.64±0.20 


8.52±0.16 


Ci OA . A 1 A 

8.89+0.14 


Nei 


8.08±0.06 


8.39+0.07 (7) 


8.42+0.06 (5) 


8.50+0.14(7) 


8.40+0.08 (9) 








8.10+0.05** 




Mgi 


7.58±0.01 


7.52±0.08(1) 


7.46+0.04(4) 


7.60±0.04(A) 




7.48+0.17 










Mgn 


7.58+0.01 


7.53±0.04 (12) 


7.43±0.04(6) 


7.40+0.05 (5) 


7. 42 ±0.02 (4) 


7.46+0.17 


7. 42 ±0.2 2 


7.59±0.22 


7.38±0.12 


7.70+0.34 


Ah 


6.49±0.01 


6.11+0.06 (2) 


6.06+0.06 (2) 
















Aln 


6.49±0.01 


6.39+0.17 (5) 


6.55+0.27 (6) 


6.38+0.38(6) 


6.14+0.08 (4) 












Aim 


6.49±0.0l 








7.00+0.38 (4) 




6.08±0.12 


6.21±0.19 


6.12±0.18 


6.31+0.15 


Sin 


7.56±0.0l 


7.58+0.19 (4) 


7.45±0.29 (3) 


7.33+0.05 (2) 


^7 ^^7 . f\ 1 1 /1\ 

7.27+0.13 (3) 


7.33+0. 17 




7. 42 ±0.2 3 


7.19±0.21 




Sim 


7.56±0.0l 








8.13+0.23 (3) 




7.26±0.19 


7. 42 ±0.2 3 


7.19±0.21 


7.50+0.15 


Pn 


5.56+0.06 








5.53+0.06 (4) 












Sn 


7.20+0.06 


7.15+0.07 '(14) 


7.07+0.08 (10) 


7. i2±0. 05 (16) 


7.05+0.09 (26) 












S in 


7.20+0.06 








7 '.08+0.04 (2) 




7.24±0.14 




6.87+0.26 




Can 


6.35+0.01 


6.31 (1) 


















TT 

OC II 


o 1 A-uA A1 
J. 1U+U.U1 


Z.J /±U. 14 


A/1 _i_n T7 


2.42 (1) 




^ 1 T-i-A 9A 

j.ij+v.zv 










Tin 


4.94+0.02 


4.89+0.13 (29) 


4.86+0.12(25) 


4.57+0.09(15) 


5.07+0.07(4) 


4.86+0.25 










Vn 


4.02+0.02 


3.57+0.06 (6) 


3.63+0.08 (6) 


3.55+0.04(3) 














Cm 


5.69+0.01 


5.62+0.08 (29) 


5.55+0.08 (24) 


5.28+0.08 (21) 


5.42+0.11(8) 


5.61+0.23 










Mnn 


5.53+0.01 


5.38+0.02(7) 


5.36+0.04(6) 


5.24(1) 


5.33(1) 


5.81+0.20 










Fei 


7.50+0.01 


7.34+0.12(21) 


7.37+0.10(13) 






7.56+0.24 










Fen 


7.50+0.01 


7.52+0.09(35) 


7.42+0.05(35) 


7.34+0.07(24) 


7.47+0.05(20) 


7.40+0.11 










Fern 


7.50+0.01 








7.45+0.07(10) 








7.36+0.20 




Nin 


6.25+0.01 


6.30+0.06 (7) 


6.17+0.04(7) 


6.00+0.01 (2) 


5.91+0.08 (4) 












Srn 


2.92+0.02 


2.37+0.04 (2) 


2.41+0.03(2) 


2.49+0.04 (2) 




2.41+0.21 










Ban 


2.22+0.02 


2.00(1) 


2.13(1) 

















a Grevesse & Sauval 1 1998 1: meteoritic abundances for non- volatile elements; b Venn ( 1995a 1995b l; c Daflon & Cunha 1 2004 1: mean values of stars from 17 clusters/associations at R ±2 kpc; 



Gummersbach et al. < 1998 »: mean values of 10 stars at /?o±2kpc; c mean from Kilian il992l [1994 1; f Rolleston et al. 1 2000 1: mean of stars from 11 clusters/associations at Rq ±2 kpc; 
* reanalysis by Venn & Przybilla <2003l - result from Venn I1995bi : 8.05±0.19\ ** Ne n 
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To conclude, effective temperatures from our approach ap- 
pear to be systematically lower than from previous model 
atmosphere analyses and higher than from previous (spec- 
tro)photometric studies. Not only is consistency now achieved 
from these indicators, but also the uncertainties in the parame- 
ter determination have been markedly reduced by bringing all 
the indicators simultaneously into agreement. This provides a 
solid basis for precise metal abundance analyses, which are dis- 
cussed next. 



6. Abundance analysis 

The chemical composition of the stars is deduced from the 
analysis of the weak line spectra of the individual elements. 
This is done by fitting non-LTE and LTE line profiles to the 
observation, as this puts tighter constraints on the modelling 
than is possible by reproducing equivalent widths alone. The 
results from the analysis of several hundred spectral lines from 
30 ionization stages of 20 chemical species are summarised in 
Appendix [A] We have also thoroughly tested the line list for 
suitability in abundance analyses of BA-stars on the main se- 
quence (Przybilla 2002 1 and recommend it for use in future 
applications. 

A comparison of non-LTE and LTE results from the in- 
dividual metal lines in 77 Leo is made in Fig.^] This super- 
giant shows the richest metal-line spectrum of all our sample 
objects. Note that a few lines from Table |A~TI without mea- 
sured equivalent widths had to be omitted for this, and results 
for a few strong 1 lines have been added from Przybilla et 
al. ^5000 1. By correctly accounting for non-LTE effects homo- 
geneous abundances are derived from all the available indica- 
tors, which show little scatter around the mean value (the de- 
rived abundance). This also implies that the microturbulent ve- 
locity is correctly chosen, and it is the same for the elements 
under study. The non-LTE analysis thus avoids the systematic 
trends of abundance with line-strength in the approximation of 
LTE, such as for N 1 or 1. In other cases, like for singly ion- 
ized sulphur, titanium and iron, the non-LTE abundances are 
systematically shifted relative to LTE. Also the statistical scat- 
ter is reduced in non-LTE, and contrary to common assump- 
tion significant non-LTE abundance corrections are found even 
in the weak line limit. These can exceed 0.3 dex for lines with 
Wa < 10 mA, in particular for the more luminous objects, where 
the non-LTE effects are amplified. The qualitative behaviour of 
non-LTE and LTE abundances is similar in all the sample stars, 
which occurs because of the reasons discussed in Sect.0] 

Mean elemental abundances for the sample stars are given 
in Table|4] along with the solar standard and results from the lit- 
erature on Galactic A-supergiants and their main sequence pro- 
genitors of early B-type at less than 2kpc distance (in order to 
avoid artefacts introduced by the Galactic abundance gradient). 
Note that we prefer solar standard abundances as summarised 
by Grevesse & Sauval ( 1998 1 over the more recent compilation 
of Asplund et al. ( 2005 1 for several reasons discussed below. 
Non-LTE results are indicated by italics. For the sample stars 
statistical uncertainties (1-cr standard deviations from the line- 
to-line scatter) are given, with the number of analysed lines in 



parenthesis. In the case of the previous studies the uncertainties 
denote 1-cr standard deviations from the star-to-star scatter. 

The present non-LTE analysis finds abundances compatible 
with the LTE study (with non-LTE corrections for a few ele- 
ments) of Venn ( 1995a 1995b 1 on less-luminous A-SGs. This 
is because the objects in that study are, on the mean, far cooler 
and of higher surface gravity than ours, which reduces the im- 
portance of non-LTE. In such a case a tailored LTE analysis 
is meaningful for most elements. Consequently, we and Venn 
recover abundances which agree within the uncertainties with 
the present-day ISM abundances (which appear to be slightly 
sub-solar), as expected. Exceptions are the carbon and nitrogen 
abundances, which reflect mixing of the surface layers of the 
supergiants by nuclear processed material, see below. 

The B main sequence stars (B-MSs) together with Hn re- 
gions define the present-day abundances in the solar neighbour- 
hood. Consequently, mutual agreement with our results exists 
for the ff-process elements with non-LTE abundances. The B- 
star abundances of helium and nitrogen are lower and the car- 
bon abundances higher than in the BA-SGs, indicating little 
to no mixing with nuclear processed matter, which is also ex- 
pected. Note that the B-MSs oxygen abundances appear to be 
slightly lower than in our BA-SGs; this may be an effect of 
small number statistics. Note also that the LTE analysis by 
Rolleston et al. (2000l finds systematically different B-MSs 
abundances than the other three B-MSs studies. Little infor- 
mation is available for the iron-group abundances, and noth- 
ing on s-process abundances in B-MSs. BA-type supergiants 
and B-MSs are therefore complementary for galactochemical 
and stellar evolution studies: the former contain information 
on the heavier element abundances and mixing patterns, while 
the latter provide stellar present-day abundances for the light 
elements. Consistency can be checked on the basis of the a- 
elements for which both star groups have features in common. 

The comparison of our sample stars with the solar standard 
(Grevesse & Sauval 1998 ) is visualised in Fig. ^] First, we 
want to discuss the results from the LTE analysis, which com- 
prises a wider variety of elements than can be done in non-LTE 
at present. For the least luminous supergiant of our sample, 
T] Leo, close to solar abundances are derived for most of the a- 
process and iron group elements. The only s-process elements 
accessible through their resonance lines, Sr n and Ba n, appear 
underabundant, as are the lightest iron group elements, scan- 
dium through vanadium. Neon appears to be overabundant and 
a discrepancy in the Al 1/11 ionization equilibrium is notewor- 
thy, while Mg 1/11 is on the verge of a discrepancy. Of the light 
elements, carbon is under- and nitrogen overabundant, with 
a discrepancy in the Ni/11 ionization balance indicated. Note 
that these LTE results are valid for the finally adopted parame- 
ters from the non-LTE analysis. A pure LTE study would indi- 
cate a slightly higher T e g from most of the ionization equilib- 
ria, implying an overall increase of the metallicity and a more 
pronounced N 1/11 discrepancy. Similar abundance patterns are 
found for the LC lab supergiant HD 111613, with several of 
the ff-processes elements (O, Ne and S) showing a tendency 
to being overabundant, as also appears to be the case for car- 
bon. For these two stars an LTE analysis appears meaningful, in 
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the sense that abundances are derived which are expected for a 
young star of the Galactic thin disk, with only a few exceptions. 

The situation changes drastically when considering an ob- 
ject of even higher luminosity, like HD 92207. Practically no 
element comes close to the solar standard, with the entire iron 
group appearing to be vastly underabundant, and the a-process 
and light elements largely overabundant. An LTE analysis sug- 
gests an abundance pattern resembling a-enhancement for this 
star, which is appropriate for the old halo population but not 
for a Population I object. This pattern is basically also derived 
for/JOri. However, at the higher atmospheric temperatures fea- 
tures from double-ionized species appear, which suggest dras- 
tic imbalances in several ionization equilibria. No parameter 
combination can be found that brings all the available ioniza- 
tion equilibria simultaneously into agreement in LTE. 

The non-LTE analysis of some of the astrophysically most 
interesting elements on the other hand reveals a rather homo- 
geneous appearance of the abundance distribution in the four 
stars. Note that we derive absolute abundances, which depend 
on physical data and not on oscillator strengths (and line broad- 
ening data) tuned to match observations, as is the case in many 
differential studies of solar-type stars. The a-process and iron 
group elements reproduce the solar abundance pattern, showing 
little scatter around slightly sub-solar mean metallicities (the 
grey bands in Fig. 1141. Among the light elements, helium and 
nitrogen are enhanced and carbon is underabundant, which is 
qualitatively well understood in the framework of the evolution 
of rotating massive stars, see the next section for details. 

The total abundance of the fusion catalysts C, N, and O 
in the sample stars correlates well with the stellar metallicity. 
Use of the solar abundances of Asplund et al. (2005) instead 
of the Grevesse & Sauval ( 1998 1 data as standard would imply 
an overabundance of CNO relative to the heavier elements in 
these young stars, in particular [O/Fe] > 0. This is inconsistent 
with other astrophysical indicators. Here, the more recent so- 
lar abundance data may imply another conflict like the recently 
unveiled solar model problem, i.e. the loss of accordance be- 
tween helioseismological measurements and the predictions of 
theoretical solar models (Bahcall et al. l2005> . 

Note that good agreement for the available non-LTE ion- 
ization equilibria is obtained, and that the statistical scatter in 
the abundance determination is typically reduced with respect 
to LTE, amounting to 0.05 to O.lOdex (1-cr standard devia- 
tions). Nitrogen and oxygen are particularly good examples. 
This comes close to the accuracy that is usually obtained in 
differential analyses of solar-type stars, however on absolute 
terms. 

The non-LTE abundance corrections vary from element to 
element and there is a well-known tendency for increased non- 
LTE corrections with luminosity (i.e. lower logg), and temper- 
ature. Typically, non-LTE abundance corrections amount to a 
factor of two to three on the mean. In some cases, the correc- 
tions can amount to factors of several tens for individual transi- 
tions, see e.g. oxygen in Fig.^] and in the most extreme cases 
they exceed a factor of a hundred. Many strong lines of Ti n 
and Fe n have been omitted in the present study, because of the 
additional computational costs and their insensitivity as abun- 
dance indicators. These also show stronger non-LTE effects as 



the corresponding weak-line spectra. This is also the case for 
the prominent Mg n ,14481 feature (Przybilla et al. l2001a> . 

The only exception from our general improvement due 
to non-LTE refinements is the Mgi/n ionization balance in 
HD 922207, which we attribute to small insufficiencies of our 
LTE model atmospheres at highest luminosity. The Mg i lines 
are highly sensitive to even the smallest changes in the mod- 
elling close to the Eddington limit, see the discussion in Sect.|3] 
In this case the Ni/ii ionization equilibrium is a better choice 
for the temperature indicator, as its behaviour is better con- 
strained - much more accurate atomic data are available for 
nitrogen than for magnesium, and the ground state ionization 
for N i is determined by the Lyman continuum instead of the 
more complicated radiation field longward of the Lyman jump, 
which is more prone to modelling uncertainties, as in the case 
of Mg i. 

The detailed non-LTE behaviour of the ionic species un- 
der study is well explained from their atomic structure, as dis- 
cussed in Sect. |4] This allows us to estimate, at least qualita- 
tively, the non-LTE behaviour of the large number of elements 
in the supergiant spectra, for which we do not have non-LTE 
spectrum synthesis at hand for the moment. Certainly, the en- 
tire singly-ionized iron group elements will show a similar non- 
LTE weakening like Tin and Fen, as they all are subject to 
non-LTE overionization. This will also be the case for Sr n and 
Ba n. A non-LTE strengthening of the Ne i lines was described 
by Sigut ( 1999) in B-MSs, which can be expected to occur to a 
similar extent also at lower T s s in the BA-SGs. This would im- 
ply a present-day neon abundance from the sample stars com- 
patible with their overall metallicity close to the solar value as 
given by Grevesse & Sauval (1998), and agreement with the 
ISM abundance as derived from the Orion nebula (Esteban et 
al. 12004b . for non-LTE abundance corrections by a factor ~2. 
Compatibility with the value of log Ne/H + 12 = 8.27 as deter- 
mined by Drake & Testa ( 2005 ) from nearby solar-like stars is 
possible for lower non-LTE abundance corrections. The solar 
value of Asplund et al. (2005) which contributes to the solar 
model problem could only be restored for unrealistically high 
non-LTE corrections by a factor ~4-5 (the Ne i lines are typ- 
ically weak). Note that an increase of the surface abundance 
of neon due to mixing with nuclear-processed matter is only 
expected in the Wolf-Rayet phase of the evolution of massive 
stars, such that BA-SGs should reflect present-day neon abun- 
dances of the ISM. 

We conclude that the derived abundance patterns from our 
LTE study are probably an artefact of the assumption of de- 
tailed equilibrium. It may be speculated that close to solar 
abundance patterns will be recovered for the four sample stars 
from a full non-LTE analysis. We will continue our efforts to 
improve the status of the modelling situation, but note that for 
many ionic species practically none of the required atomic data 
needed to construct reliable non-LTE model atoms are avail- 
able at present. 

Another important aspect of the non-LTE abundance anal- 
ysis are constraints for the systematic uncertainties. This in- 
cludes tests for the response of the spectrum synthesis to mod- 
ifications of the stellar parameters as well as to variations 
of the atomic data that is used for assembling the non-LTE 



N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



23 



model atoms. Such tests have been done for Ci/ii, Ni/n, Oi 
and Mgi/n by Przybilla et al. ( EUOUl 12001 al b) and Przybilla & 
Butler (2001 1. When we scale these results for the stellar pa- 
rameter uncertainties of the present work and consider contin- 
uum placement in our high-S/N spectra to be a negligible prob- 
lem we can conclude that the systematic uncertainties of our 
non-LTE abundance determination are of the order ~0.10dex 
(1-cr-error). This is from quadratic summation over all uncer- 
tainty contributors, assuming them to be independent. As this 
is not strictly the case for the stellar parameters (e.g. T^/logg 
have to be varied simultaneously to retain consistency of a 
spectroscopic indicator with observation, see e.g. Fig. I10i. the 
true systematic uncertainties may be slightly smaller. The given 
uncertainty has consequently to be viewed as a conservative es- 
timate. Note that the main contribution to the systematic uncer- 
tainty of abundance determinations in BA-SGs comes from the 
uncertainties in the stellar parameters. The typical uncertainties 
in oscillator strengths, photoionization and collisional excita- 
tion cross-sections from quantummechanical ab-initio compu- 
tations are of the order 10-20% and therefore contribute only 
little to the total error budget. Thoroughly constraining the sys- 
tematic uncertainties is an extremely computationally intensive 
task, requiring hundreds of model runs with modified parame- 
ters, the problem thus being ~2 orders of magnitude more com- 
plex than the abundance determination itself. We therefore as- 
sume the systematic abundance uncertainties for S n/m, Ti n and 
Fen also to amount to ~0.10dex (1-cr-error). The atomic data 
used for assembling these model atoms is of similar quality 
(perhaps slightly lower because the atoms are more complex) 
than in the lighter elements. The response to stellar parameter 
changes can also be expected to be similar, as deduced from a 
few exploratory test calculations. 

We omit a detailed comparison of our non-LTE abundance 
determination with previous work from the literature. For C i/ii, 
N i/n, O i and Mg i/ii this has been done elsewhere (Przybilla et 
al. lMTQllMJTal b: Przybilla & Butler l2UoTV The main conclu- 
sions are that in general the results from the other studies are 
close to our findings, or origins of the differences have been 
found, e.g. different stellar parameters, see Sect. 15. 41 Major dis- 
crepancies were found only for N i, where the previous stud- 
ies gave systematically lower nitrogen abundances, by a fac- 
tor ~ 2. This has been traced to use of approximate collision 
data in the previous non-LTE model atoms, but see also Venn 
& Przybilla (2003 1. Contrary to previous studies, our non-LTE 
analysis for these elements yields homogeneous results from 
all available spectral lines. Non-LTE studies of S n/m, Ti n and 
Fe ii in BA-SGs have not been reported in the literature so far. 

Finally, a comparison of the high-resolution observations 
of the luminous Galactic supergiant HD 92207 with our spec- 
trum synthesis is made in Fig.^] The model can successfully 
reproduce the line spectrum with high accuracy over extended 
wavelength regions, except for a few features which are strong 
enough to be affected by sphericity effects and the stellar wind. 
An obvious example for this is Hy. Some complications in lo- 
calised spectral regions occur because of interstellar absorption 
bands, which are not modelled. Finally, CCD defects and in- 
sufficiently corrected cosmic ray hits can render some regions 
useless for the analysis, e.g. near Ti n /14501 in the present case. 



Table 5. Surface abundance data relevant for stellar evolution 



Element 


?/Leo 


HD 111613 


HD 92207 


/JOri 


Y 


0.37+0.04 


0.32+0.04 


0.35+0.04 


0.38+0.04 


N/C 


2.93+1.38 


1.89+0.69 


0.97+0.09: 


2.68+0.52 


N/O 


0.37+0.10 


0.44+0.12 


0.25+0.05 


0.45+0.09 


[CNO/H] 


-0.06+0.05 


-0.08+0.04 


-0.04+0.05: 


±0.00+0.04 


[M/H] 


-0.04+0.03 


-0.11+0.03 


-0.09+0.07 


-0.06+0.10 



In high-resolution data these are easily recognised. However at 
intermediate resolution, in particular at low S/N, special care 
has to be taken to identify the artefacts and to exclude them 
from the analysis. 

7. Evolutionary status 

Massive stars provide the dominant contribution to the momen- 
tum and energy budget of the interstellar matter. They are the 
main drivers for the dynamical and chemical evolution (as sites 
of nucleosynthesis) of the interstellar environment, and conse- 
quently for the evolution of galaxies. Understanding the evo- 
lution of massive stars is a prerequisite for such further stud- 
ies. Massive star evolution is an active field of research (see 
e.g. Maeder & Meynet (2000 1 for an overview) where the over- 
all problem has been solved in many respects but where many 
details are still not conclusively explained because of the un- 
certainty in several free parameters of the input physics. These 
have to be constrained by observation. In particular quantifying 
the role of stellar winds and rotation as a function of metallic- 
ity requires further efforts, as well as transport phenomena as 
manifested in surface abundance anomalies, and the existence 
and extent of blue loops. BA-type supergiants are of special 
interest in this regard, as they allow these problems to be inves- 
tigated beyond the Galaxy and its nearest satellites, in a much 
wider range of different galactic environments. In the following 
we wish to discuss our four sample stars in the context of their 
evolutionary status, putting special emphasis on areas where 
our analysis technique allows for an improvement on previous 
work. 

First, reduced uncertainties in the stellar parameters (see 
Table |2) allow the objects to be located precisely in the 
Hertzsprung-Russell diagram (HRD, Fig. I16i. Here, the major 
remaining uncertainties arise from the constraints of an indirect 
distance determination. Improvements on this can only be ex- 
pected from future astrometric satellite missions. The objects 
cover a wide range in luminosity within the supergiant regime, 
and include two of the most luminous and massive BA-SGs 
analysed so far. Evolutionary and spectroscopic masses are in 
excellent agreement, indicating good consistency between the 
stellar properties inferred from the present analysis and those 
predicted by theory. In particular the integrated mass-loss his- 
tory of these (near-)solar metallicity objects is apparently well 
described. 

Then, additional information on the evolutionary status of 
the stars comes from abundance patterns of the light elements. 
The reaction rates for the CNO-cycle result in an accumulation 
of nitrogen at the cost of carbon and (to lesser extent) oxygen, 
with helium being the fusion product. If transport mechanisms 
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are active throughout the stellar envelope, these fusion prod- 
ucts may be mixed into the stellar atmosphere. The observed 
abundance patterns are summarised in Table [5] (by mass frac- 
tion for Y, N/C and N/O). The most sensitive indicator, the N/C 
ratio, is also discussed in the HRD diagram in Fig. ^] All 
objects show enrichments of helium and nitrogen, depletion 
of carbon and practically solar oxygen. The sum of the CNO 
abundances [CNO/H] correlates well with the stellar metallic - 
ity [M/H], consistent with the catalyst role of these elements in 
the main fusion cycle on the main sequence. The comparison 
with the predictions of evolution computations by Meynet & 
Maeder (2003) shows good qualitative agreement, while dis- 
crepancies remain in the details. In particular, helium enrich- 
ment traces the increasing N/C ratio well, while oxygen deple- 
tion obviously does not occur to the predicted extent. 

In this scenario the three more luminous supergiants appear 
to have evolved directly from the main sequence, at initial ro- 
tational velocities larger than 300 kms 1 (ySOri, HD 111613) 
and at a lower rate, respectively, in the case of HD 92207. 
The high helium abundance and N/C ratio (the highest of the 
sample) of 77 Leo on the other hand are indicative for the first 
dredge-up and thus a blue loop scenario for this object. Note 
that this classification can only be made for the individual ob- 
jects with statistical significance because of the reduced abun- 
dance uncertainties. However, additional efforts are required to 
further improve the significance of the N/C ratios, which are 
currently still impaired by the remaining uncertainties in the 
carbon abundance determination. 

8. Analyses at intermediate resolution 

In order to analyse intermediate-resolution data {R 1 000- 
5 000) whole spectral regions have to be modelled at once, as 
the individual spectral features are typically no longer resolved 
by the instrument. The spectrum synthesis technique imple- 
mented here allows for such analyses, as practically all spectral 
features in BA-SGs of significant strength (more than a few 
mA) are included. In particular the spectral region between the 
Balmer jump and ~5000A - the classical region for analyses 
of hot stars - is well covered. Thus, an estimate of the over- 
all stellar metallicity and abundances of individual elements 
from several unblended lines of intermediate strength can be 
obtained from a comparison of the observed with synthesised 
spectra. The feasibility of this approach is tested in the follow- 
ing. Note that spectral regions containing lines of intermediate 
strength (Wa < 300 mA) are best suited for an analysis, as these 
are strong enough to produce a noticeable signal even at mod- 
est S/N, but are also weak enough to be of photospheric origin. 
In BA-SGs several ionic species give rise to such lines, typi- 
cally N 1, 1, Mg 11, Si 11, Ti n, Cr 11, Fe 11 and in the hotter stars 
He 1 as well. The three iron group elements hereby dominate 
the line spectrum, whereas He 1, N 1 and the three a-elements 
contribute only a few distinct features. 

A test for the applicability of our spectrum synthesis at in- 
termediate resolution is made in the following, based on the 
most luminous object of our sample, HD 92207. The transition 
from the high resolution comparison to intermediate resolution 
is made in Fig.^J where the observed and the synthetic spec- 



trum from Fig. have been artifically degraded to 5 A reso- 
lution. Such data can be expected e.g. from observations with 
FORS1 (FOcal Reducer and low-dispersion Spectrograph) at 
the VLT. The excellent agreement between model and obser- 
vation is preserved, except for a few spectral features already 
discussed above (Sect.|6j. However, observational data at such 
a high S/N will typically not be available. In order to con- 
sider a more realistic case, the observed spectrum is further 
degraded to a S/N ^ 50 in Fig. ED Additionally, test calcula- 
tions for scaled abundances and modified effective temperature 
are displayed. While the surface gravity can still be determined 
at high accuracy from the higher Balmer lines, testing for the 
temperature sensitivity of the synthetic spectra becomes nec- 
essary as ionization equilibria are unavailable for a direct Jeff- 
determination (typically, one of the ionic species shows only 
weak lines). Instead, empirical (metallicity-dependent) spectral 
type-r e ff calibrations have to be applied (see e.g. Kudritzki 
et al. 120031 Evans & Howarth 2003 ), increasing the uncer- 
tainties, or a spectrophotometric approach, which requires fur- 
ther investigation (see Sect. 15.31 . The tests show that within 
a considerable r e fj-interval the systematic effects on the over- 
all spectrum synthesis are rather small. Modifications of the 
abundances on the other hand show a much larger effect, and 
we conclude that our method allows the stellar metallicity to 
be constrained within approximately +0.2 dex (1-cr uncertain- 
ties). This is considerable larger than the statistical uncertain- 
ties from the high-resolution analysis but sufficient for mean- 
ingful applications in extragalactic stellar astronomy. In fact, 
such uncertainties are characteristic for the status of contem- 
porary abundance work on BA-SGs in the literature. From a 
few features, like Fen /L1435 1,4508 and the pure Fen blends 
around 4520 and 4580 A, even the iron abundance can be de- 
rived in this wavelength region. In a similar manner, one can 
use the regions around He 1 AA4Q26 and 447 1 to determine the 
helium abundance. More chemical species become accessible 
in other wavelength regions. 

Note that the microturbulent velocity can also not be di- 
rectly derived from intermediate resolution spectra, as the 
whole range from the weak lines to those of intermediate 
strength is not available for the analysis. However, our high 
resolution analyses and data from the literature indicate that 
microturbulence values of £=8+2kms~ 1 are appropriate for 
BA-SGs at luminosities log L/Lq> 5.0. We assume that this 
empirical finding holds for analogous objects in the galaxies 
beyond the Local Group, for which only the next generation of 
extremely large telescopes will allow for high-resolution spec- 
troscopy and direct microturbulence determinations. Within the 
typical uncertainties in no marked systematic effects for the 
metallicity determination are found. 

To conclude, the present spectrum synthesis technique in 
combination with sophisticated non-LTE model atoms allows 
for reliable metallicity determinations for BA-type supergiants 
on the basis of intermediate-resolution spectra. With suitable 
spectra available, abundances for several elements of gen- 
eral astrophysical relevance can also be obtained. This opens 
up BA-SGs as versatile tools for extragalactic stellar astron- 
omy. First applications to metallicity analyses of A-type su- 
pergiants well beyond the Local Group, in the field spiral 
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galaxy NGC 3621 and the Sculptor group spiral NGC 300, have 
already been discussed by Przybilla ( 2002 1 and Bresolin et 

al. ( BOOllliOOa . 

9. Summary and conclusions 

We have presented a hybrid non-LTE spectrum synthesis tech- 
nique for the quantitative analysis of BA-type supergiants and 
related objects of lower luminosity and tested it on high-quality 
observations of four bright Galactic objects. The statistical 
equilibrium and radiative transfer problem for several of the 
astrophysically most important elements is solved on the ba- 
sis of classical line-blanketed LTE atmospheres. These have 
been shown to be sufficient to obtain consistent results from the 
available spectroscopic indicators. A prerequisite for this is that 
secondary atmospheric parameters such as helium abundance, 
metallicity and microturbulence are accounted for consistently, 
because of the high sensitivity of the model atmospheres and 
line formation to the details of the computations close to the 
Eddington limit. It has been shown that the stellar parameters of 
BA-SGs can be determined with unprecedented precision: the 
internal accuracy of the methods allows T e g to be constrained 
to better than 1-2% and logg to 0.05-0.10dex. Non-LTE com- 
putations with sophisticated model atoms reduce random errors 
and remove systematic trends in abundance analyses. The l-cr 
errors in the determination of absolute abundances amount to 
typically 0.05-0. lOdex and ~0.10dex for statistical and sys- 
tematic uncertainties, respectively. This comes close to the ac- 
curacy that is usually achieved in differential studies of solar- 
type stars in the contemporary literature. Contrary to common 
assumption, significant non-LTE abundance corrections by a 
factor 2-3 can occur even for the weakest lines in BA-SGs, and 
considerably larger corrections for the stronger lines. Because 
of this, a purely LTE analysis of highly luminous supergiants 
fails to obtain meaningful results, while the less-luminous su- 
pergiants can in principle be addressed in such a way, if re- 
duced accuracy can be tolerated. LTE analyses thereby tend to 
systematically overestimate the abundances of the light and a- 
process elements, and to underestimate the abundances of the 
iron group species. The non-LTE analysis implies close to solar 
abundance for the heavier elements in the sample objects, and 
patterns characteristic for mixing with nuclear-processed mat- 
ter in the case of the light elements. These indicate a blue-loop 
scenario for one of the sample stars because of first dredge-up 
abundance ratios, while the other objects were restricted to the 
blue part of the HRD. Using the full potential of our spectrum 
synthesis, we have shown that almost the entire visual and near- 
IR spectra of BA-SGs can be accurately reproduced, implying 
the the technique is also suitable for analyses of intermediate- 
resolution spectra. Thus BA-type supergiants have become a 
diagnostic tool for quantitative stellar spectroscopy beyond the 
Local Group. 

Appendix A: Spectral line analysis 

In this appendix we provide details on our spectral line analy- 
sis, as a basis for further applications. Table |A~TI summarises 
our line data and the results from the abundance analysis 



for our sample stars. It is available in electronic form only. 
The first columns give the line wavelength A (in A), excita- 
tion energy of the lower level x ( m eV), adopted oscillator 
strength log gf, an accuracy indicator and the source of the gf 
value. Most of the g/-data is retrieved from the atomic spec- 
tra database (V2.0) of the National Institute of Standards and 
Technology (http://physics.nist.gov/cgi-bin/AtData/mainjsd I. 
Then, for each object the measured equivalent width W,\ 
(in mA) is tabulated, followed by the derived abundance 
log s = log (X/H) +12. In cases with an entry for the non-LTE 
abundance correction A logs = log £nlte - log £lte this de- 
notes the non-LTE abundance, else the LTE abundance. The 
equivalent widths have been measured by direct integration 
over the spectral lines, deviating from the usual approach of 
fitting the observed profiles by a Gaussian. For high-S/N ob- 
servations of BA-SGs as in the present study this helps to avoid 
systematic uncertainties, as the line profiles are potentially sub- 
ject to asymmetries imposed by the velocity field at the base 
of the stellar wind. In some cases metal lines are situated in 
the wings of the Balmer lines. Then their equivalent width is 
measured with respect to the local continuum, indicated by 
'S(Wa)'. Abundances are determined from a best match of the 
spectrum synthesis to the observed line profiles, and not the 
equivalent widths, see the discussion in Sect. [6] Equivalent 
widths are used as auxiliary diagnostics for the determination 
of microturbulent velocity (see e.g. Fig.1131. Spectrum synthe- 
sis also allows blended features to be used for abundance deter- 
minations, where a W,\ measurement is hampered. These cases 
are marked by a sole 'S' in Table IA.fi For Hei only non-LTE 
abundances are derived because of the potential of this major 
atmospheric constituent to change atmospheric structure and 
thus the stellar parameter determination, see Sect. 13.21 Entries 
in italics are disregarded when computing abundance averages 
(Table|Hl. 
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20.96 


-2.28 


B 


W JU 


















^SA7 


20.96 


-2.75 


B 


W iu 


















QQ7 1 7Q 

jo / i./y 


21.22 


-1.85 


C 


u/cr 
W ju 


S 


11.24 










132 


11.13 


QQOA CLA 

jyZO. j4 


21.22 


-1.65 


A 




S 


11.11 








Al 11 17 

4Z 11.1/ ... 


S 


11.17 


'iQfi.A I'X 
OV04. / 


20.62 


-1.30 


A 


W iu 














S(152) 


11.24 


zinna 9 A 


21.22 


-1.47 


C 


W iU 












zl7 11 1 A 
4 / 11. 14 


156 


11.16 


4UZJ.V0 


21.22 


-2.58 


B 


W JU 


S 


11.17 










S 


11.19 


4UZ0. lo 


20.96 


-2.63 


A 


Wou 


S 


11.14 




c 1 1 (19 
11 .uz 




1 oa 11 m 

1ZO 11. 1U ... 


S 


11.16 


/inoA 1 o 
4uzo. ly 


20.96 


-0.63 


A 


W iu 


















4UZ0.ZU 


20.96 


-0.85 


A 


W JU 




















20.96 


-1.32 


A 


W 1 V 1 


















/l 1 on 8 1 
41ZU.01 


20.96 


-1.74 


B 


W iU 


S 


11.17 




z / 1 1 .u / 




11 17 
Jy 11.1/ ... 


118 


11.20 


41ZU.0Z 


20.96 


-1.96 


B 


W iU 


















/i 1 on qq 


20.96 


-2.44 


B 


u/cr 
W iu 


















A\A^X 7 A 

414j. 1 


21.22 


-1.20 


B 


W JU 














S 


11.23 


/I 1 Q7 


21.22 


-2.34 


A 




11 


11.19 




in 11 1 a 

1U 1 1. LH 






43 


11.17 


4jo 1 .yj 


21.22 


-0.88 


A 


W iu 


S 


11.17 




K 11 11 
1 1.1 j 




SO 111C 

oU 1 1. lo ... 


S 


11.14 


44j /. DJ 


21.22 


-2.03 


B 


W JU 


10 


11.15 










56 


11.19 


AA1 1 /I7 
44 / 1.4 / 


20.96 


-0.20 


A 




S 


11.16 




s 11.00 




C 1 1 HA 

i> 1 l.Uu 


S 


11.23 


4471.49 


20.96 


-0.42 


A 


WSG 


















4471.68 


20.96 


-0.90 


A 


WSG 


















4713.14 


20.96 


-1.23 


B 


WSG 


43 


11.18 




40 11.12 




58 11.17 


164 


11.15 


4713.16 


20.96 


-1.45 


B 


WSG 


















4713.38 


20.96 


-1.93 


B 


WSG 


















4921.93 


21.22 


-0.44 


A 


WSG 


70 


11.19 




c 1 1 no 

S 11.08 




S 11.19 






5875.60 


20.96 


-1.52 


A 


WSG 


161 


11.24 




142 11.05 




a 11 nn 

248 1 1.09 






5875.61 


20.96 


0.48 


A 


WSG 


















5875.63 


20.96 


-0.34 


A 


WSG 


















5875.64 


20.96 


0.14 


A 


WSG 


















5875.97 


20.96 


-0.22 


A 


WSG 


















Ci: 

7113.18 


8.65 


-0.77 


B- 


WFD 


5 


8.03 


-0.21 












9078.29 


7.48 


-0.58 


B 


WFD 


S 


7.88 


-0.32 


S 8.05 


-0.48 








9088.51 


7.48 


-0.43 


B 


WFD 


S 


8.01 


-0.52 












9111.81 


7.49 


-0.30 


B 


WFD 


S 


7.82 


-0.40 


S 8.15 


-0.62 








Cn: 

4267.00 


18.05 


0.56 


C+ 


WFD 


36 


8.17 


-0.18 


40 8.18 


-0.28 


66 8.33 -0.32 


158 


8.10 -0.41 


4267.26 


18.05 


0.74 


C+ 


WFD 


















5132.95 


20.70 


-0.21 


B 


WFD 














9 


8.20 -0.19 


5133.28 


20.70 


-0.18 


B 


WFD 


















6151.27 


20.84 


-0.15 


D 


F 














S 


8.16 +0.01 


6151.53 


20.84 


0.02 


D 


F 


















6578.05 


14.45 


-0.03 


B 


WFD 


S(12) 


8.12 


-0.15 


28 8.25 


-0.28 








6582.88 


14.45 


-0.33 


B 


WFD 


S 


8.00 


-0.13 


17 8.30 


-0.26 









Ni: 



3830.43 


10.69 


-1.39 


C+ 


WFD 


S 


8.37 


-0.11 




















5999.43 


11.60 


-1.41 


c+ 


WFD 


5 


8.49 


-0.13 




















6008.47 


11.60 


-1.11 


c+ 


WFD 


10 


8.52 


-0.17 




















7423.64 


10.33 


-0.71 


B+ 


WFD 


76 


8.40 


-0.48 


87 


8.45 


-0.54 


49 


8.24 


-0.55 


32 


8.49 


-0.51 


7442.30 


10.33 


-0.38 


B+ 


WFD 


130 


8.39 


-0.71 


152 


8.41 


-0.66 


80 


8.19 


-0.61 


59 


8.50 


-0.60 


7468.31 


10.34 


-0.19 


B+ 


WFD 


165 


8.34 


-0.98 


197 


8.46 


-0.81 


125 


8.26 


-0.69 


90 


8.45 


-0.68 


7898.98 


12.36 


0.02 


C 


WFD 


30 


8.38 


-0.31 


19 


8.35 


-0.32 














7899.28 


12.36 


-0.91 


C 


WFD 



























Table Al. (cont.) 



N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



rjheo HD 111613 HD92207 /JOri 

/1(A) x (eV) log gf Acc. Src. VK,i(mA) logs Aloge H^(mA) logs A logs W A (mA) logs Aloge M^(mA) log£ Aloge 
Nii 



8567.74 


10.68 


-0.66 


B 


WFD 


S(95) 


8.43 


-0.44 


S(85) 


8.41 


-0.47 


S(40) 


8.18 


-0.52 








8594.00 


10.68 


-0.33 


B 


WFD 


S(84) 


8.36 


-0.93 


S 


8.32 


-0.69 














8629.24 


10.69 


0.08 


B 


WFD 


S(184) 


8.14 


-0.86 


S(209) 


8.18 


-0.69 














8655.88 


10.69 


-0.63 


B 


WFD 


S(72) 


8.29 


-0.51 


S(57) 


8.32 


-0.46 














8680.28 


10.34 


0.35 


B+ 


WFD 


S(373) 


8.48 


-1.71 


S(532) 


8.53 


-1.57 


S 


8.29 


-1.27 


S 


8.62 


-1.38 


8683.40 


10.33 


0.09 


B+ 


WFD 


S(312) 


8.51 


-1.40 


S(408) 


8.51 


-1.22 


S 


8.26 


-0.96 


S 


8.57 


-0.95 


8686.15 


10.33 


-0.31 


B+ 


WFD 




















s 


8.60 


-0.67 


8703.25 


10.33 


-0.32 


B+ 


WFD 


211 


8.55 


-0.94 


239 


8.49 


-0.63 


123 


8.27 


-0.64 


s 


8.46 


-0.58 


8711.70 


10.33 


-0.23 


B+ 


WFD 


228 


8.47 


-1.02 


266 


8.54 


-0.77 


155 


8.21 


-0.68 


s 


8.51 


-0.61 


8718.84 


10.34 


-0.34 


B+ 


WFD 


198 


8.50 


-0.92 


233 


8.47 


-0.62 


122 


8.29 


-0.63 


s 


8.45 


-0.56 


8728.90 


10.33 


-1.07 


B+ 


WFD 


S(63) 


8.40 


-0.34 


S(63) 


8.39 


-0.34 














9028.92 


11.60 


-0.13 


B 


WFD 


s 


8.39 


-0.27 


s 


8.26 


-0.36 














9045.88 


12.36 


0.44 


B 


WFD 


112 


8.38 


-0.49 


S 


8.38 


-0.48 


44 


8.25 


-0.56 


s 


8.41 


-0.61 


9049.49 


12.36 


-0.86 


B 


WFD 


82 


8.38 


-0.43 








S 


8.31 


-0.47 


s 


8.48 


-0.53 


9049.89 


12.36 


0.28 


B 


WFD 


























Nii: 
3955.85 


18.47 


-0.81 


B 


WFD 




















17 


8.50 


-0.18 


3995.00 


18.50 


0.21 


B 


WFD 


10 


8.32 


-0.15 


10 


8.36 


-0.18 


15 


8.28 


-0.29 


65 


8.42 


-0.45 


4447.03 


20.41 


0.23 


B 


WFD 




















27 


8.52 


-0.26 


4601.48 


18.46 


-0.43 


B+ 


WFD 




















26 


8.43 


-0.24 


4607.15 


18.46 


-0.51 


B+ 


WFD 




















23 


8.48 


-0.22 


4613.87 


18.46 


-0.67 


B+ 


WFD 




















19 


8.39 


-0.21 


4630.54 


18.48 


0.09 


B+ 


WFD 




















S 


8.49 


-0.41 


4643.09 


18.48 


-0.36 


B+ 


WFD 




















S 


8.54 


-0.30 


4788.14 


20.65 


-0.36 


B 


WFD 




















8 


8.53 


-0.09 


4803.29 


20.67 


-0.11 


B 


WFD 




















15 


8.61 


-0.12 


5045.10 


18.46 


-0.41 


B+ 


WFD 




















25 


8.57 


-0.35 


5666.63 


18.47 


-0.05 


A 


WFD 




















30 


8.49 


-0.37 


5676.02 


18.46 


-0.37 


A 


WFD 




















19 


8.50 


-0.29 


5679.56 


18.48 


0.25 


A 


WFD 




















51 


8.60 


-0.54 


5686.21 


18.47 


-0.55 


A 


WFD 




















15 


8.55 


-0.29 


5710.77 


18.48 


-0.52 


A 


WFD 




















16 


8.59 


-0.32 


Oi: 

3947.29 


9.15 


-2.10 


B 


WFD 


35 


8.68 


-0.04 


34 


8.68 


-0.11 


34 


8.82 


-0.18 


13 


8.75 


-0.22 


3947.48 


9.15 


-2.24 


B 


WFD 


























3947.59 


9.15 


-2.47 


B 


WFD 


























4368.19 


9.52 


-2.67 


B- 


WFD 


34 


8.75 


-0.11 














1 1 


8.83 


-0.22 


4368.24 


9.52 


-1.96 


B- 


WFD 


























4368.26 


9.52 


-2.19 


B- 


WFD 


























4654.12 


10.74 


-2.16 


C+ 


WFD 


12 


8.83 


-0.05 




















4654.56 


10.74 


-1.93 


C+ 


WFD 


























4772.45 


10.74 


-1.92 


c+ 


WFD 


14 


8.80 


-0.04 


16 


8.72 


-0.08 














4772.91 


10.74 


-1.70 


c+ 


WFD 


























4773.75 


10.74 


-1.55 


c+ 


WFD 


10 


8.68 


-0.03 




















4967.38 


10.74 


-1.63 


c+ 


WFD 


47 


8.78 


-0.07 


41 


8.76 


-0.10 


34 


8.81 


-0.21 








4967.88 


10.74 


-1.41 


c+ 


WFD 


























4968.79 


10.74 


-1.26 


c+ 


WFD 


























5329.10 


10.74 


-1.24 


c+ 


WFD 


50 


8.78 


-0.11 


47 


8.63 


-0.22 


67 


8.78 


-0.23 


32 


8.75 


-0.25 


5329.68 


10.74 


-1.02 


c+ 


WFD 


























5330.73 


10.74 


-0.87 


c+ 


WFD 


39 


8.78 


-0.11 


34 


8.69 


-0.22 














6046.23 


10.99 


-1.76 


c+ 


WFD 


S 


8.83 


-0.09 




















6046.44 


10.99 


-1.54 


c+ 


WFD 


























6046.49 


10.99 


-2.24 


c+ 


WFD 


























6155.96 


10.74 


-1.36 


B+ 


WFD 


120 


8.80 


-0.14 


110 


8.70 


-0.17 


183 


8.73 


-0.33 


113 


8.73 


-0.35 


6155.97 


10.74 


-1.01 


B+ 


WFD 


























6155.99 


10.74 


-1.12 


B+ 


WFD 


























6156.74 


10.74 


-1.49 


B+ 


WFD 



























30 

Table Al. (cont.) 



N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



z/Leo HD 111613 HD92207 /30ri 

/1(A) x (eV) log gf Acc. Src. ff,i(mA) loge A loge W A (m A) logs A loge W A (mA) loge A logs IV^(mA) loge A logs 



Oi: 



6156.76 


10.74 


-0.90 


B+ 


WFD 
























6156.78 


10.74 


-0.69 


B+ 


WFD 
























6158.15 


10.74 


-1.84 


B+ 


WFD 


89 


8.80 


-0.17 


89 8.66 


-0.20 














6158.17 


10.74 


-1.00 


B+ 


WFD 
























6158.19 


10.74 


-0.41 


B+ 


WFD 
























6453.60 


10.74 


-1.29 


C+ 


WFD 


38 


8.77 


-0.09 


S 8.75 


-0.19 


S 


8.88 


-0.30 


S 


8.78 


-0.30 


6454.44 


10.74 


-1.07 


C+ 


WFD 
























7001.90 


10.99 


-1.49 


B 


WFD 


S 


8.80 


-0.15 


S 8.68 


-0.21 


S 


8.69 


-0.35 


S 


8.83 


-0.32 


7001.92 


10.99 


-1.01 


B 


WFD 
























7002.17 


10.99 


-2.66 


B 


WFD 
























7002.20 


10.99 


-1.49 


B 


WFD 
























7002.23 


10.99 


-0.74 


B 


WFD 
























7002.25 


10.99 


-1.36 


B 


WFD 
























On 

3954.36 


23.42 


-0.40 


B 


WFD 


















9 


8.85 


-0.24 


4349.43 


23.00 


0.06 


B 


WFD 


















S(12) 


8.84 


-0.33 


4366.89 


23.00 


-0.35 


B 


WFD 


















7 


8.85 


-0.27 


4641.81 


22.98 


0.05 


B 


WFD 


















15 


8.84 


-0.27 


4661.63 


22.98 


-0.28 


B 


WFD 


















6 


8.79 


-0.23 


Nei: 
































5852.49 


16.85 


-0.49 


B 


S 


5 


8.38 








10 


8.40 




23 


8.37 




6030.00 


16.67 


-1.04 


B- 


WSG 


















7 


8.31 




6074.30 


16.67 


-0.50 


B 


S 


7 


8.30 




11 8.50 




10 


8.31 




25 


8.28 




6096.16 


16.67 


-0.31 


B 


s 


1 1 


8.33 




11 8.43 




17 


8.45 




39 


8.41 




6143.06 


16.62 


-0.10 


B 


S 


13 


8.51 




13 8.41 




S 


8.47 




59 


8.52 




6163.59 


16,72 


-0.62 


B 


s 


















22 


8.33 




6266.50 


16.72 


-0.37 


B 


s 


8 


8.37 




9 8.34 




15 


8.49 




32 


8.45 




6334.43 


16.62 


-0.32 


B 


s 


















38 


8.46 




6383.00 


16.67 


-0.23 


B 


s 


S 


8.42 




















6402.20 


16.62 


0.33 


B 


s 


19 


8.40 




24 8.40 




52 


8.71 










6506.50 


16.62 


-0.03 


B 


s 












22 


8.65 




55 


8.49 




Mgi: 
































3829.36 


2.71 


-0.21 


B 


WSM 


S 


7.49 


+0.06 


S 7.47 


+0.13 


S 


7.62 


+0.19 








3832.30 


2.71 


0.27 


B 


WSM 


S 


7.41 


+0.19 


S 7.44 


+0.18 


S 


7.60 


+0.23 








3838.29 


2.72 


0.49 


B 


WSM 


S 


7.44 


+0.14 


S 7.42 


+0.17 


S 


7.63 


+0.23 








4702.99 


4.35 


-0.42 


C+ 


BMZ 


12 


7.64 


+0.14 


















5172.68 


2.71 


-0.38 


B 


WSM 


S 


7.53 


+0.15 


















5183.60 


2.72 


-0.16 


B 


WSM 


90 


7.51 


+0.12 


102 7.50 


+0.13 


23 


7.54 


+0.14 








5528.41 


4.35 


-0.40 


C+ 


BMZ 


6 


7.60 


+0.17 


















Mgn: 
































3848.21 


8.86 


-1.49 


c 


T 


S(59) 


7.52 


+0.05 


S(46) 7.42 


+0.02 


31 


7.41 


+0.11 


14 


7.45 


+0.06 


3848.34 


8.86 


-2.40 


D 


T 
























3850.39 


8.86 


-1.74 


C 


T 


S 


7.54 


-0.02 


S 7.48 


+0.01 


S 


7.48 


+0.06 








4384.64 


10.00 


-0.79 


c+ 


WSM 


S 


7.54 


-0.02 


















4390.51 


10.00 


-1.71 


D 


WSM 


95 


7.54 


-0.02 


103 7.47 


±0.00 


52 


7.36 


+0.05 


37 


7.41 


+0.03 


4390.57 


10.00 


-0.53 


D 


WSM 
























4427.99 


10.00 


-1.20 


C+ 


WSM 


32 


7.56 


+0.02 


















4433.99 


10.00 


-0.90 


c+ 


WSM 


49 


7.55 


+0.02 


















4739.59 


11.57 


-0.66 


c+ 


T 


33 


7.50 


+0.01 


21 7.39 


-0.01 


9 


7.37 


+0.04 


12 


7.42 


+0.02 


4739.71 


11.57 


-0.77 


c+ 


T 
























4851.08 


11.63 


-0.42 


c 


CA 


S(24) 


7.51 


+0.02 


S(17) 7.45 


-0.01 


S(10) 7.40 


-0.03 


S(7) 


7.41 


+0.02 


5401.54 


11.63 


-0.08 


c 


CA 


S 


7.57 


-0.03 


















6545.97 


11.63 


0.41 


c 


CA 


S(75) 


7.61 


-0.11 


S(73) 7.38 


-0.12 














7877.05 


10.00 


0.39 


c+ 


WSM 


S 


7.45 


-0.43 


















7896.04 


10.00 


-0.30 


c+ 


WSM 


S 


7.48 


-0.49 


















7896.37 


10.00 


0.65 


c+ 


WSM 

























N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



Table Al. (cont.) 



A (A) 


^(eV) logg/ Acc. 


Src. 


77 Leo 




HD 111613 


HD 92207 


/30ri 




W^(mA) loge 


A loge 


W A (mA) loge A loge 




A log s 


W^(mA) loge 


A loge 


Ah: 


























3944.OI 


0.00 


-0.64 


C+ 


WSM 


34 6.15 




35 6.10 












3961.52 


0.01 


-0.34 


c+ 


WSM 


S^43 , l 6 06 




S 6.01 












Aln: 


























4663.05 


10.60 


-0.29 


c 


WSM 


54 6.41 




S 6.26 


57 6.16 




38 


6.06 




5593.30 


13.26 


0.41 


D 


WSM 


9 6.09 
















6226.19 


13.07 


0.03 


D- 


WSM 


8 6.45 
















6731 6? 


13.07 


-0.09 


D- 


WSM 






20 6.39 


16 6.07 










6231.75 


13.07 


0.39 


D 


WSM 




















13.08 


-1.27 


E 


WSM 


77 ft 40 




7fi ft ^7 


21 6.12 










6243 70 


13.08 


-0.09 


D- 


WSM 


















62A331 


13.08 


0.65 


D 


WSM 


















7042.08 


11.32 


0.34 


C+ 


WSM 






S 6.89 


110 6.96 




63 


6.19 




7056.71 


11.32 


0.11 


C+ 


WSM 






S 6.86 


S 6.75 




S 


6.22 




7471.41 


13.65 


0.74 


D 


WSM 












11 


6.07 




8640.71 


11.82 


-0.02 


C 


WSM 


S(29) 6.50 




43 6.58 


S 6.22 










Al m: 


























4^17 ^7 


17.81 


0.42 


c+ 


WSM 












18 


6.67 




4528.95 


17.82 


-0.28 


c+ 


WSM 












25 


6.69 




4529.19 


17.82 


0.67 


c+ 


WSM 


















5696 60 


15.64 


0.24 


B 


WSM 












57 


7.39 




5722.73 


15.64 


-0.07 


B 


WSM 












36 


7.26 




Sin: 


























3853 66 


6.86 


— 1.0U 


1 

b 


WSM 












124 


TIC 

f.LJ 




407^ 4S 


9.84 


-1.29 


X 


MEL 


SHM 7 41 




S 7 ^4 


S 7.29 










4621.42 


12.53 


-0.54 


D 


WSM 


27 7.43 




S 7.23 












4621.70 


12.53 


-1.68 


D 


WSM 


















4621.72 


12.53 


-0.39 


D 


WSM 


















5669.56 


14.20 


0.32 


X 


MEL 








S 7.36 










5957.56 


10.07 


-0.35 


D 


WSM 


S 7.69 




S 7.78 






S 


7.26 




5978.93 


1U.U / 


U.UD 




WSM 


127 7.80 










1 1 7 


7 41 




Sim: 


























4552.62 


19.02 


0.29 


C+ 


WSM 












69 


8.37 




4567.84 


19.02 


0.07 


C+ 


WSM 












45 


8.10 




4574.76 


19.02 


-0.41 


c+ 


WSM 












23 


7.91 




Pn: 

S7S3 48 


11.02 


0.32 


D 


WSM 












S 


5.47 




5296 08 


10.80 


-0.16 


D 


WSM 












7 


5.55 




5425 88 

JT^Z/J . OO 


10.80 


0.18 


D 


WSM 












17 


5.60 




6043 08 

\j\J'-T*J .\JO 


10.80 


0.42 


D 


WSM 












19 


5.48 




Sn: 

4153.07 


15.90 


0.62 


D- 


WSM 


16 7.18 


-0.10 


18 7.11 -0.17 


23 7.20 


-0.15 


40 


7.09 


-0.13 


4162.67 


15.94 


0.78 


D- 


WSM 


17 7.09 


-0.11 


S 6.99 -0.15 


S 7.12 


-0.08 


50 


6.96 


-0.10 


4269.73 


16.09 


-0.12 


E 


WSM 












11 


7.12 


-0.07 


4278.51 


16.09 


-0.11 


E 


WSM 












11 


7.09 


-0.07 


4282.59 


16.10 


-0.01 


E 


WSM 












13 


7.08 


-0.06 


4294.40 


16.13 


0.58 


D- 


WSM 












32 


6.90 


-0.18 


4552.41 


15.07 


-0.13 


D 


WSM 








17 7.25 


-0.33 








4716.27 


13.62 


-0.41 


D 


WSM 


7 7.15 


-0.06 


13 7.10 -0.11 


18 7.11 


-0.17 


37 


7.03 


-0.15 


4815.55 


13.67 


0.09 


D 


WSM 


17 7.08 


-0.08 


21 7.00 -0.13 


38 7.12 


-0.31 


77 


7.04 


-0.27 


4885.65 


14.00 


-0.61 


D 


WSM 


6 7.29 


-0.10 








21 


7.07 


-0.15 


5009.57 


13.62 


-0.28 


D 


WSM 


11 7.21 


-0.10 




27 7.18 


-0.23 


47 


7.11 


-0.18 


5014.04 


14.07 


0.10 


D 


WSM 


14 7.20 


-0.10 




S 7.13 


-0.29 


S 


7.19 


-0.46 


5032.43 


13.67 


0.27 


D 


WSM 












94 


6.91 


-0.53 


5212.27 


15.07 


-0.62 


D 


WSM 








S 7.07 


-0.29 


38 


6.92 


-0.16 


5212.62 


15.07 


0.32 


D 


WSM 
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z/Leo HD 111613 HD92207 /30ri 

/1(A) x (eV) log gf Acc. Src. ff,i(mA) loge A loge H^(mA) logs A loge W A (mA) loge A logs W^fmA) loge A logs 



Sn: 


































5320.72 


15.07 


0.50 


D 


WSM 


10 


7.12 


-0.12 


13 


7.09 


-0.12 


19 


7.00 


-0.42 


57 


6.89 


-0.41 


5345.71 


15.07 


0.36 


D 


WSM 




















48 


6.93 


-0.32 


5428.66 


13.58 


-0.13 


D 


WSM 


13 


7.22 


-0.08 


14 


7.13 


-0.16 








62 


7.11 


-0.36 


5432.80 


13.62 


0.26 


D 


WSM 




















108 


7.18 


-0.51 


5453.86 


13.67 


0.48 


D 


WSM 


25 


7.10 


-0.18 


34 


7.07 


-0.27 


63 


7.17 


-0.55 








5473.61 


13.58 


-0.18 


D 


WSM 




















56 


7.11 


-0.32 


5509.71 


13.62 


-0.14 


D 


WSM 


13 


7.20 


-0.11 














59 


7.13 


-0.34 


5564.96 


13.67 


-0.32 


D 


WSM 








10 


7.18 


-0.17 


17 


7.08 


-0.32 


40 


7.06 


-0.33 


jdUd.1 j 
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L) 
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1 

10 


7 (Yl 

/.U3 
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2X) 


A GO 


n 1 q 

— u.iy 


33 


n no 
/.UZ 


— U.33 








5616.63 


13.66 


-0.64 


D 


WSM 














9 


7.06 


-0.22 


20 


7.03 


-0.14 


5647.02 


14.00 


0.04 


D 


WSM 


12 


7.19 


-0.13 








32 


7.25 


-0.39 








SAAO no 
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VV Olvl 










7 OA 


-fl 91 
— U.Z 1 


zz 


7 1 9 
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U.Z 1 


^1 
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— n 19 


SAA4. 77 


1 ^ AA 


—0 9S 

U.Z, J 


LJ 


VV OlVl 


u 


7 in 

/ . 1U 


— n 1 9 

U. 1 Z 








1 1 
1 1 


7 fn 


— n 99 

— U.ZZ 


^8 

JO 


U.7U 


— n 9 1 

u.z 1 


A^1 1 AQ 


IzL 11 


U.J 1 


c 


OP 




















90 
Zu 


7 1 n 


in 


6397.36 


14.16 


-0.41 


C 


OP 




















28 


7.11 


-0.17 


6398.01 


14.15 


-0.81 


c 


OP 


























6413.71 


14.16 


-0.81 


c 


OP 




















9 


7.17 


-0.18 


s in- 

JyZo.jO 


1 8 ^9 


— U. 1 / 


p 






















<; 
J 


7 1 1 


n iq 

U. I7 


tiJJ.JU 


1 8 9zl 


u.JU 


LJ 


VV OlVl 




















19 
1 z 


7 


-f) 9zt 


Can: 
5019.97 


7.51 


-0.28 


D 


WSM 


21 


6.31 






















Sen: 


































4246.82 


0.32 


0.28 


D 


KB 


36 


2.41 




72 


2.45 




19 


2.42 










4415.56 


0.60 


-0.64 


X 


KB 


S 


2.64 






















5526.79 


1.77 


0.06 


D 


KB 


8 


2.67 




15 


2.83 
















Tin: 
3900.56 


1.13 


-0.45 


D 


MFW 


140 


4.95 


+0.30 


S(225) 


4.82 


+0.24 


96 


4.85 


+0.41 








o r\ 1 o /i o 

3913.48 


1.12 


-0.53 


D 


MFW 


132 


4.97 


+0.13 


202 


4.80 


+0.25 


85 


4.81 


+0.34 


9 


5.05 


+0.30 


4028.36 


1.89 


-1.00 


D 


MFW 


32 


4.95 


+0.33 


58 


5.03 


+0.23 


S 


4.91 


+0.37 








4163.63 


2.59 


-0.40 


D 


MFW 


55 


5.07 


+0.28 


91 


5.03 


+0.27 


S 


4.89 


+0.28 








4171.92 


2.60 


-0.56 


D 


MFW 


54 


5.11 


+0.28 




















4287.88 


1.08 


-2.02 


D- 


MFW 


15 


5.02 


+0.26 




















4290.22 


1.16 


-1.12 


D- 


MFW 


68 


4.92 


+0.28 








35 


4.94 


+0.35 








4290.35 


2.06 


-1.53 


X 


KB 


























4294.09 


1.08 


-1.11 


D- 


MFW 








124 


4.89 


+0.22 


46 


4.92 


+0.27 








4300.06 


1.18 


-0.77 


D- 


MFW 


1 17 


5.13 


+0.33 








76 


5.01 


+0.33 








4301.92 


1.16 


-1.16 


D- 


MFW 


42 


4.81 


+0.32 


73 


4.90 


+0.27 














4312.87 


1.18 


-1.16 


D- 


MFW 








91 


4.93 


+0.24 


29 


4.93 


+0.34 








4314.97 


1.16 


-1.13 


D- 


MFW 


59 


4.88 


+0.29 


S 


4.87 


+0.25 














4330.24 


2.04 


-1.51 


D 


MFW 


S(7) 


4.84 


+0.28 


S(30) 


4.92 


+0.21 














4330.72 


1.18 


-2.04 


D- 


MFW 


S(7) 


4.90 


+0.29 




















4394.02 


1.22 


-1.59 


D- 


MFW 


15 


4.65 


+0.27 


S 


4.75 


+0.24 














4395.00 


1.08 


-0.66 


D- 


MFW 


114 


4.86 


+0.38 


S 


4.80 


+0.18 


82 


4.90 


+0.44 








4399.79 


1.24 


-1.27 


D- 


MFW 


39 


4.88 


+0.28 


S 


4.92 


+0.26 


20 


4.88 


+0.40 








4417.72 


1.16 


-1.43 


D- 


MFW 


38 


4.96 


+0.27 




















4443.78 


1.08 


-0.70 


D- 


MFW 


91 


4.77 


+0.37 


S 


4.73 


+0.24 


61 


4.76 


+0.41 


5 


4.95 


+0.16 


4450.50 


1.08 


-1.45 


D- 


MFW 


30 


4.71 


+0.25 


52 


4.74 


+0.23 














4468.52 


1.13 


-0.60 


D- 


MFW 


93 


4.64 


+0.30 


154 


4.65 


+0.30 


46 


4.67 


+0.44 


S 


4.96 


+0.27 


4501.27 


1.12 


-0.75 


D- 


MFW 


79 


4.72 


+0.33 


147 


4.67 


+0.28 














4563.77 


1.22 


-0.96 


D- 


MFW 


78 


4.95 


+0.28 


132 


4.90 


+0.24 


45 


4.91 


+0.29 








4571.96 


1.57 


-0.53 


D- 


MFW 


111 


4.94 


+0.24 


176 


4.90 


+0.20 


68 


4.90 


+0.26 


8 


5.09 


+0.32 


4779.98 


2.05 


-1.37 


D- 


MFW 


17 


4.98 


+0.23 


38 


5.05 


+0.25 














4798.53 


1.08 


-2.43 


D 


MFW 


5 


4.92 


+0.27 


7 


4.81 


+0.17 














4805.09 


2.06 


-1.10 


D- 


MFW 


39 


4.95 


+0.27 


55 


4.98 


+0.19 
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Fig. 14. Results from the elemental abundance analysis (relative to the solar composition) for our sample stars. Filled symbols 
denote non-LTE, open symbols LTE results. The symbol size codes the number of spectral lines analysed - small: 1 to 5, medium: 
6 to 10, large: more than 10. Boxes: neutral, circles: single-ionized, diamonds: double-ionized species. The error bars represent 
lcr-uncertainties from the line-to-line scatter. The grey shaded area marks the deduced metallicity of the objects within lcr-errors. 
Non-LTE abundance analyses reveal a striking similarity to the solar abundance distribution, except for the light elements which 
have been affected by mixing with nuclear-processed matter. 
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Fig. 15. Comparison of our spectrum synthesis (thin red line) with the high-resolution spectrum of HD 92207 (full line). The 
major spectral features are identified, short vertical marks designate Fe n lines. The location of the diffuse interstellar band 
around 4430 A is indicated. With few exceptions, excellent agreement between theory and observation is found (Hy is affected 
by the stellar wind). 
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Fig. 17. Same as Fig.^J however artifically degraded to 5 A resolution. The excellent agreement between theory and observation 
is preserved. Deviations are found only for Hy, Mg n /1448 1 (see text) and in the regions of the interstellar band around 4430 A 
and near Ti n ^4501 (due to a CCD defect, cf. also Fig.fBl 
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Fig. 18. Same as Fig.^J but artifically degraded to a S/N ^ 50. The comparison with synthetic spectra for scaled abundances and 
modified temperatures (thin red lines, as indicated) shows that the models react sensitively to changes of the metal abundances 
only. Consequently, stellar parameters estimated from the spectral classification (T e g) and from fitting of the higher Balmer lines 
(log g) suffice to constrain the stellar metallicity to ~0.2 dex. The lower set of spectra in each panel has been shifted by -0.2 units. 
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z/Leo HD 111613 HD92207 /30ri 

/1(A) ^ (eV) log g/ Acc. Src. W A (mA) logs A logs W A (mA) loge A logs W A (mA) loge A logs IV^(mA) logs A logs 



Tin: 






















4911.18 


3.12 


-0.34 


D 


MFW 


21 


4.60 +0.20 


34 


4.64 +0.15 






5129.16 


1.89 


-1.39 


D- 


MFW 


20 


4.96 +0.33 


31 


4.95 +0.22 






5188.68 


1.58 


-1.21 


D- 


MFW 


32 


4.91 +0.30 


50 


4.93 +0.24 


12 


4.75 


5336.78 


1.58 


-1.70 


D- 


MFW 


15 


4.90 +0.29 


S 


4.93 +0.26 






Vn: 

3916.42 


1.43 


-1.06 


B 


MFW 


9 


3.61 


22 


3.73 






3951.97 


1.48 


-0.74 


B 


MFW 


17 


3.51 


31 


3.55 


S 


3.51 


3997.13 


1.48 


-1.20 


C 


MFW 






17 


3.69 






4005.71 


1.82 


-0.46 


D 


MFW 


27 


3.58 


S 


3.65 


14 


3.55 


4023.39 


1.80 


-0.52 


X 


KB 


18 


3.51 


29 


3.54 






4035.63 


1.79 


-0.62 


X 


KB 


16 


3.56 


27 


3.62 


S 


3.58 


4183.44 


2.05 


-0.95 


X 


KB 


8 


3.65 
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3979.51 


5.65 


-0.73 


X 


KB 


S(40) 


5.61 


53 


5.57 


28 


5.40 


4037.97 


6.49 


-0.56 


X 


KB 


22 


5.56 


22 


5.49 






4070.84 


6.49 


-0.75 


X 


KB 


27 


5.72 


22 


5.61 






4086.13 


3.71 


-2.42 


X 


KB 


S(ll) 


5.55 


13 


5.51 






4242.36 


3.87 


-1.33 


X 


KB 


96 


5.79 


138 


5.67 


63 


5.35 


4252.63 


3.84 


-2.02 


X 


KB 










20 


5.32 


4261.91 


3.87 


-1.53 


X 


KB 


73 


5.61 


99 


5.60 


38 


5.30 


4269.28 


3.85 


-2.17 


X 


KB 


27 


5.69 


39 


5.66 


17 


5.29 


4275.57 


3.86 


-1.71 


X 


KB 


54 


5.68 


76 


5.65 


31 


5.31 


4284.19 


3.86 


-1.86 


X 


KB 


45 


5.70 


64 


5.64 


28 


5.34 


4555.01 


4.07 


-1.38 


D 


MFW 


65 


5.60 










4558.66 


4.07 


-0.66 


D 


MFW 


159 


5.75 






134 


5.36 


4588.22 


4.07 


-0.63 


D 


MFW 


146 


5.54 






111 


5.14 


4592.07 


4.07 


-1.22 


D 


MFW 


68 


5.46 


95 


5.42 


42 


5.20 


4616.64 


4.07 


-1.29 


D 


MFW 


60 


5.52 


84 


5.45 


30 


5.17 


4618.82 


4.07 


-1.11 


D 


MFW 


111 


5.77 






78 


5.39 


4634.10 


4.07 


-1.24 


D 


MFW 


92 


5.72 


138 


5.70 






A O 1 1 ~> ,1 

4812.34 


3.86 


-1.99 


X 


KP 


30 


5.59 


38 


5.52 


15 


5.22 


4824.12 


3.87 


-0.96 


X 


KP 


127 


5.61 






95 


5.20 


4836.22 


3.86 


-1.94 


X 


KP 


30 


5.58 


45 


5.57 


16 


5.27 


4848.24 


3 86 


-1.14 


x 


KB 




J.J7 . . . 




S S6 


79 

/ y 




4876.41 


3.86 


-1.46 


D 


MFW 


S(88) 


5.57 


S(119) 


5.50 


56 


5.24 


4884.58 


3.86 


-2.08 


D 


MFW 


S(21) 


5.60 


31 


5.57 






4901.62 


6.49 


-0.83 


X 


KB 


12 


5.54 










5237.33 


4.06 


-1.16 


D 


MFW 


83 


5.51 


127 


5.49 






5246.76 


3.71 


-2.45 


D 


MFW 


13 


5.56 


19 


5.59 






5274.96 


4.05 


-1.29 


X 


KB 






S 


5.48 


S 


5.37 


5308.46 


4.07 


-1.85 


D 


KB 


27 


5.57 


32 


5.46 


15 


5.28 


5313.61 


4.07 


-1.65 


D 


MFW 






65 


5.57 


24 


5.31 


5334.87 


4.05 


-1.56 


X 


KB 






53 


5.43 


25 


5.15 


5420.91 


3.76 


-2.36 


D 


MFW 


13 


5.59 


22 


5.51 






5502.07 


4.17 


-1.99 


D 


MFW 


20 


5.61 










5508.63 


4.16 


-2.11 


D 


MFW 


18 


5.66 










Mnn 






















4206.37 


5.40 


-1.57 


X 


KB 


11 


5.39 


10 


5.35 






4251.72 


6.18 


-1.06 


X 


KB 


11 


5.40 


S 


5.31 






4326.64 


5.40 


-1.25 


X 


KB 


S(16) 


5.36 


s 


5.33 






4478.64 


6.65 


-0.95 


X 


KB 


8 


5.39 


11 


5.34 






4755.73 


5.40 


-1.24 


X 


KB 


17 


5.39 


31 


5.41 






5297.00 


9.86 


-0.21 


X 


KB 


8 


5.33 


7 


5.39 


S 


5.24 


5297.03 


9.86 


0.43 


X 


KB 














5297.06 


9.86 


0.62 


X 


KB 















7 5.47 



33 5.45 

25 5.26 

4 5.24 

15 5.51 



21 5.51 

S(14) 5.42 
S(9) 5.48 



6 5.33 



Table Al. (cont.) 



N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



z/Leo HD 111613 HD92207 /30ri 

/1(A) x (eV) log gf Acc. Src. ff,i(mA) loge A loge H^(mA) logs A logs W A (mA) loge A logs IV^(mA) loge A logs 



Fei: 



3815.84 


1.48 


0.30 


C+ 


FMW 


S(48) 


7.19 




73 


7.24 
















3859.91 


0.00 


-0.71 


B+ 


FMW 


51 


7.26 




80 


7.28 
















3922.91 


0.05 


-1.65 


B+ 


FMW 


7 


7.35 




13 


7.44 
















3927.92 


0.11 


-1.59 


C 


FMW 


14 


7.28 




16 


7.45 
















3948.77 


3.27 


-0.34 


C 


FMW 


7 


7.45 






















4045.81 


1.48 


0.28 


B+ 


FMW 


60 


7.24 




79 


7.29 
















4063.59 


1.56 


0.07 


C+ 


FMW 


42 


7.23 




53 


7.26 
















4071.74 


1.61 


-0.02 


B+ 


FMW 


35 


7.22 




S 


7.28 
















4118.54 


3.57 


0.28 


C 


FMW 


S(8) 


7.34 




S(H) 


7.57 
















4132.06 


1.61 


-0.65 


D- 


FMW 


14 


7.36 






















4181.75 


2.83 


-0.18 


D- 


FMW 


7 


7.42 




s 


7.49 
















4202.03 


1.48 


-0.71 


B+ 


FMW 


16 


7.43 






















4260.47 


2.40 


-0.02 


D 


FMW 


24 


7.51 




24 


7.46 
















4271.15 


2.45 


-0.35 


B+ 


FMW 


13 


7.48 




S 


7.42 
















4271.76 


1.48 


-0.16 


B+ 


FMW 


37 


7.33 






















4278.23 


3.37 


-1.74 


C 


FMW 


43 


7.40 




52 


7.33 
















4383.54 


1.48 


0.20 


B+ 


FMW 


59 


7.47 






















4404.75 


1.56 


-0.14 


B+ 


FMW 


35 


7.26 




40 


7.34 
















4415.12 


1.61 


-0.62 


B+ 


FMW 


15 


7.36 






















5615.64 


3.33 


-0.14 


C+ 


FMW 


7 


7.55 






















6862.48 


4.55 


-1.57 


c 


FMW 


11 


7.10 






















Fen: 
3779.58 


2.53 


-3.78 


X 


KB 








47 


7.31 


+0.06 


30 


7.30 


+0.16 








3781.51 


4.48 


-2.77 


X 


KB 








44 


7.34 


+0.06 


22 


7.28 


+0.15 








3783.35 


2.27 


-3.16 


X 


KB 








171 


7.38 


+0.07 


87 


7.20 


+0.20 


S 


7.60 


+0.25 


3922.00 


9.13 


-1.07 


X 


KB 


8 


7.41 


±0.00 


7 


7.37 


+0.09 














3938.29 


1.67 


-3.89 


D 


FMW 


77 


7.47 


+0.09 


S 


7.43 


+0.07 














3945.21 


1.70 


-4.25 


D 


FMW 


47 


7.44 


+0.06 


68 


7.43 


+0.06 


37 


7.43 


+0.25 








4031.44 


4.73 


-3.12 


X 


KB 


12 


7.46 


+0.07 


15 


7.33 


+0.04 














4041.64 


5.57 


-3.13 


X 


KB 


7 


7.48 


+0.01 


30 


7.40 


+0.06 














4122.64 


2.58 


-3.38 


D 


FMW 


S(87) 


7.64 


+0.11 


S(130) 


7.53 


+0.05 


62 


7.41 


+0.24 








4124.79 


2.54 


-4.20 


D 


FMW 


S(23) 


7.55 


+0.09 


S(30) 


7.54 


+0.04 


11 


7.39 


+0.18 








4233.17 


2.58 


-2.00 


C 


FMW 




















117 


7.44 


+0.40 


4273.32 


2.70 


-3.34 


D 


FMW 


87 


7.52 


+0.10 


114 


7.37 


+0.04 


54 


7.30 


+0.29 


11 


7.51 


+0.40 


4296.57 


2.70 


-3.01 


D 


FMW 


127 


7.55 


+0.12 


182 


7.38 


+0.04 


103 


7.27 


+0.24 


21 


7.49 


+0.37 


4303.17 


2.70 


-2.49 


C 


FMW 




















47 


7.36 


+0.41 


a a orv i r\ 

4489.19 


2.83 


-2.97 


D 


FMW 


119 


7.54 


+0.10 


170 


7.44 


+0.04 


86 


7.25 


+0.23 


17 


7.47 


+0.42 


a a m An 

4491.40 


2.86 


-2.70 


C 


FMW 


140 


7.43 


±0.00 


205 


7.32 


+0.04 








28 


7.42 


+0.42 


4508.28 


2.86 


-2.31 


D 


KB 




















63 


7.40 


+0.40 


4515.34 


2.84 


-2.48 


D 


FMW 




















44 


7.39 


+0.41 


4520.23 


2.81 


-2.60 


D 


FMW 




















42 


7.47 


+0.40 


4522.63 


2.84 


-2.11 


C 


KB 




















86 


7.46 


+0.40 


4541.52 


2.86 


-3.05 


D 


FMW 


110 


7.59 


+0.05 


167 


7.45 


+0.06 


88 


7.30 


+0.18 


21 


7.55 


+0.41 


4555.89 


2.83 


-2.32 


D 


KB 




















56 


7.40 


+0.39 


4576.33 


2.84 


-3.04 


D 


FMW 


116 


7.49 


+0.01 


163 


7.47 


+0.05 


87 


7.30 


+0.20 


20 


7.52 


+0.40 


4580.06 


2.58 


-3.73 


X 


KB 


39 


7.37 


+0.08 


S 


7.31 


+0.04 














4582.84 


2.84 


-3.10 


C 


FMW 


89 


7.35 


+0.04 


S 


7.37 


+0.06 


S 


7.41 


+0.22 








4620.51 


2.83 


-3.28 


D 


FMW 


76 


7.40 


+0.04 


102 


7.36 


+0.06 


50 


7.22 


+0.20 








4629.34 


2.81 


-2.37 


D 


FMW 




















S 


7.36 


+0.42 


4635.33 


5.96 


-1.65 


D- 


FMW 


86 


7.70 


+0.01 




















4656.97 


2.89 


-3.63 


E 


FMW 


53 


7.52 


+0.06 


68 


7.44 


+0.06 


40 


7.42 


+0.19 








4666.75 


2.83 


-3.33 


D 


FMW 


74 


7.45 


+0.03 


108 


7.43 


+0.06 


56 


7.33 


+0.23 


9 


7.51 


+0.42 


4670.17 


2.58 


-4.10 


E 


FMW 


34 


7.58 


+0.09 


47 


7.54 


+0.06 














4731.44 


2.89 


-3.36 


D 


FMW 


86 


7.78 


+0.08 


117 


7.55 


+0.05 


49 


7.40 


+0.20 


10 


7.53 


+0.33 


4993.35 


2.81 


-3.65 


E 


FMW 


49 


7.52 


+0.06 


65 


7.35 


+0.05 














5278.94 


5.91 


-2.41 


X 


KB 


13 


7.48 


+0.03 


S 


7.40 


±0.00 














5325.56 


3.22 


-3.22 


X 


KB 


69 


7.54 


+0.06 


91 


7.41 


+0.04 


45 


7.32 


+0.21 


7 


7.51 


+0.29 



40 

Table Al. (cont.) 



N. Przybilla et al.: Quantitative Spectroscopy of BA-type Supergiants 



77 Leo 



HD 111613 



HD 92207 



ySOri 



A (A) 


^r(eV) log gf Acc. 


Src. 


W^(mA) logs 


Aloge 


H^(mA) logs 


Aloge 


W,i(mA) loge 


Aloge 


W A (mA) loge Aloge 


Fen: 
































5362.87 


3.20 


-2.74 


X 


KB 














114 


7.35 


+0.19 






5425.25 


3.20 


-3.36 


D 


FMW 


60 


7.58 


+0.03 


87 


7.45 


+0.05 


38 


7.34 


+0.22 






5427.83 


6.72 


-1.66 


X 


KB 


29 


7.51 


-0.01 


















5534.83 


3.24 


-2.93 


D 


FMW 


106 


7.56 


-0.01 


159 


7.52 


+0.02 


79 


7.37 


+0.22 


15 


7.45 +0.38 


6147.74 


3.89 


-2.72 


X 


KB 


75 


7.52 


+0.01 


101 


7.45 


+0.03 


63 


7.41 


+0.21 






6149.26 


3.89 


-2.72 


X 


KB 


77 


7.53 


±0.00 


102 


7.38 


+0.03 


64 


7.40 


+0.20 






6175.15 


6.22 


-1.98 


X 


KB 


35 


7.55 


+0.04 


39 


7.49 


+0.05 












6179.38 


5.57 


-2.60 


X 


KB 


19 


7.52 


+0.02 


21 


7.38 


+0.03 












6238.39 


3.89 


-2.63 


X 


KB 


86 


7.50 


+0.02 


110 


7.40 


+0.01 


65 


7.35 


+0.22 


19 


7.64 +0.33 


6416.92 


3.89 


-2.85 


D 


FMW 


72 


7.61 


+0.04 


80 


7.47 


+0.04 


48 


7.41 


+0.17 






6432.68 


2.89 


-3.74 


D 


FWM 


43 


7.55 


+0.05 


48 


7.42 


+0.03 












6433.81 


6.22 


-2.47 


X 


KB 


8 


7.44 


+0.09 


















Fe m: 
































4005.04 


11.58 


-1.76 


X 


KB 




















12 


7.42 


4371.13 


8.24 


-2.99 


X 


KB 




















15 


7.49 


4382.51 


8.25 


-3.02 


X 


KB 




















13 


7.45 


4395.76 


8.26 


-2.60 


X 


KB 




















28 


7.48 


4419.60 


8.24 


-2.22 


X 


KB 




















38 


7.32 


4431.02 


8.25 


-2.57 


X 


KB 




















27 


7.46 


5063.42 


8.65 


-2.95 


X 


KB 




















10 


7.53 


5086.70 


8.66 


-2.59 


X 


KB 




















S 


7.51 


5127.39 


8.66 


-2.22 


X 


KB 




















50 


7.53 




8.66 


-2.56 


X 


JVD 
























5156.11 


8.64 


-2.02 


X 


KB 




















36 


7.34 


Nin: 


































4.03 


-1.88 


X 


1/13 


92 


6.37 




S 


6.13 










S 


5.86 


4015.47 


4.03 


-2.42 


X 


KB 


46 


6.30 




56 


6.16 










11 


5.99 


4067 03 


4.03 


-1.84 


X 


KB 


95 


6.20 




134 


6.09 










35 


5.83 


4187.85 


4.03 


-2.68 


X 


KB 








42 


6.21 




30 


5.99 




5 


5.96 


4192.07 


4.03 


-3.06 


X 


KB 


15 


6.29 




17 


6.20 




S 


6.00 








4244.78 


4.03 


-3.11 


X 


KB 


13 


6.28 




S 


6.17 














4279.22 


6.62 


-2.00 


X 


KB 


10 


6.36 




















4362.10 


4.03 


-2.72 


X 


KB 


29 


6.31 




s 


6.20 




S 


6.00 








Srn: 
































4077.71 


0.00 


0.15 


X 


FW 


42 


2.34 




64 


2.39 




s 


2.52 








4215.52 


0.00 


-0.17 


X 


FW 


24 


2.40 




37 


2.43 




s 


2.46 








Ban: 
































4554.03 


0.00 


0.14 


X 


D 


5 


2.00 




9 


2.13 















accuracy indicators - uncertainties within: AA: 1%; A: 3%; B: 10%; C: 25%; D: 50%; E: larger than 50%; X: unknown 

sources of g/-values - BMZ: Butler et al. 1 1993 1; CA: Coulomb approximation (Bates & Damgaard 1949 1; D: Davidson et al. 1 1992 1; F: Fernley 
et al. (available from Topbase); FMW: Fuhr et al. (1988 1; KB: Kurucz & Bell 119951 : MEL: Mendoza et al. (available from Topbase); 
MFW: Martin et al . IT9881 : S: Sigut (19991 : T: Taylor (available from Topbase); W FD: W iese et al. fT99o1 : WSG: Wiese et al. il966l ': 
WSM: Wiese et al. ( 1969 1*; when available'*', improved gf- values from Fuhr & Wiese 1 1998 1 are favoured 

sources for Stark broadening parameters - Hi: Stehle & Hutcheon 119991 . Vidal et al. [1973 1; Hei: Barnard et al. 1 1969 1974 1, Dimitrijevic & 
Sahal-Brechot (l990l ;Ci: Griem (19741 , Cowley fiWTV Cn: Griem {T%4\ [19741 . Cowley (WTH ; Ni/w. Griem (T954lfT974l . Cowley 119711 : 
Oi/ii: Cowley (197T1 ; Nei: Griem (19741 , Cowley (19711 : Mgi: Dimitrij evic & Sahal-Brechot (19961 Cowley flWII : Mgn: Griem IT9S41 
f!974l, Co wley iTWII; Ali: Grie m 119741 . Cowley (19711; Aln: Griem n964lfl974l . Cowley IT57TI; Aim: Cowle y fl97H : Sin: Lanzet 
al. IT9881 . Griem Il974l . Cowley JTWII ; Si in: Cowley (197T1 ; P 11: Cowley ITWII ; S n/m: Cowley (197T1 ; Ca 11: Griem (19741 . Cowley fTWII : 
Sc - Ni: Cowley fPml : Srn: Cowley (TWTI : Ban: Dimitrijevic & Sahal-Brechot (19971 



